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ABSTRACT 



The Korean Air Force has determined that repairables management is one of the 
areas in which attention could be expected to lead to substantial improvement in the 
efficient management of defense resources and in maintaining an adequate level of force 
effectiveness. This thesis reviews various inventory models for the management of re- 
pairable items and develops an inventory’ model for the Korean Air Force. It discusses 
the characteristics of each model, and, identifies and explains the differences in each 
model with respect to assumptions, objectives, constraints, and optimization methods. 
Also this research evaluates the proposed model using sample data sets. 
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I. INTRODUCTION 



A. BACKGROUND 

The objective of the Korean Air Force supply system is to ensure the desired level 
of operational availability for its population of repairable equipment and maintain suf- 
ficient stock levels for replacement components and repair parts to support the mainte- 
nance of existing weapon systems and equipment. Also, the Korean Air Force supply 
system is responsible for stocking a number of items which are failure prone and have 
been designated as repairable. Because of engineering and economic considerations, at- 
tempts are made by Korean Air Force repair activities to restore these items to service- 
able condition whenever they fail. 

The repairable inventorv' system differs from that of the consumable inventory sys- 
tem in two ways. First, the repairable system contains two distinct inventories, one of 
which has items that are in a Ready-For-lssue (RFI) state and another which has items 
that are in a Non-Ready-For-Issue (NRFI) state. The first inventor\' contains items 
which are usable and the second contains items that must be repaired before they can 
be used. Second, the RFI inventorx’ is made up of a mixture of new items and items that 
have been used, failed, repaired, and are ready to be used again. This is one reason why 
repairable inventory management is a complicated process. [Ref 1: p. 2] 

The Korean Air Force Logistics Command (KAFLC) tries to determine an appro- 
priate stock level for repairables based on a primitive mathematical model. This model 
estimates stock levels using rules of thumb based on expert opinion. However, it is ver}’ 
difficult to determine an optimal stock level for a repairable without using a more so- 
phisticated inventory management model. The result of the lack of such a model is that 
the Korean Air Force supply system tends to maintain relatively high stock levels for 
repairable components. Although, repairable items account for only a small percentage 
of the quantity of items stocked throughout the Korean Air Force, they account for a 
large portion of dollars invested in inventory. The characteristics of high cost and high 
essentiality that these items typically possess make efficient inventoiy control difficult 
and extremely critical. 

Historically, most inventoiy’ models have been developed for the private sector 
where the profit motive is important. Such models consider the various average annual 
variable costs associated with managing inventories and strive to minimize the sum of 
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these costs. The typical relevant costs in most inventor}' models are ordering cost, 
holding cost and stockout cost. Certainly, these cost parameters are ver\’ important in 
the common concept of inventory management. [Ref. 2; p. 1] 

However, in militar}' supply systems, it is hard to estimate the costs associated with 
stockouts. In addition, the military supply system is not interested in profit 
maximization or cost minimization. Instead, such an organization is usually interested 
in maximizing the ability of its forces to respond to any threat. An objective which 
maximizes some measure of readiness with given resources is therefore appropriate. 

B. PURPOSE OF RESEARCH 

The effective accomplishment of each flying mission is dependent upon having suf- 
ficient numbers of aircraft ready to fly and to perform at their fullest capability. To 
support this goal, the Korean Air Force emphasizes an extensive system of supply and 
maintenance capabilities. The Korean Air Force also tries to maximize the availability 
of aircraft by quickly identifying malfunctioning parts, removing them, and rapidly in- 
stalling replacements which have been positioned at the base supply department. 

The purpose of this thesis is the development of an inventor}' model for the man- 
agement of repairable items at the wholesale level by the Korean Air Force Logistics 
Command (KAFLC). This will be accomplished by the analysis of various repairable 
inventory management models which have been developed and applied to the manage- 
ment of inventories of repairable items in the United States militar}'. And also, a nu- 
merical analysis will be accomplished based on a proposed model. Finally the inventory 
system parameters, and management concept will be studied. 

C. PREVIEW 

Chapter II of this thesis presents an overview of the current Korean Air Force Lo- 
gistics Command (KAFLC) repairable items management process. This overview will 
explain the organization of the Korean Air Force Logistics Command, describe the 
present general methodology used for repairable items management and the current 
overall approach to the problem. 

Chapter III gives a detailed description of the specific mathematical inventory 
models in use today for the management of repairables. Some explanation of their 
mathematical approach in solving the inventory problem is given. 

In Chapter IV, the development of a model for the Korean Air Force repairables 
inventory management system is described. Model formulation and solution techniques 
are provided. 
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Chapter V will show numerical results from a computer program for comparing the 
proposed model and describe major components based on the findings in Chapters 111 
and IV. 

Finally. Chapter VI contains a summar>’ and recommendations for future study 
based on the proposed repairables inventor}' management model. 
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II. REPAIRABLES AND THE KOREAN AIR FORCE LOGISTICS 

SYSTEM 



A. INTRODUCTION 

The Korean Air Force didn't use any kind of repairables inventory’ management 
model until the late 1960's. In the beginning of the 1970's, the first inventor)’ manage- 
ment concept was introduced by the U.S Air Force. The economic inventor)’ manage- 
ment model was used during this period. 

The model adopted in the 1970's has been revised systematically several times since 
its adoption by the Korean Air Force. However, the usefulness of this inventory model 
has decreased as the weapon systems used in the Korean Air Force have become more 
sophisticated, and the overall size of the Korean Air Force has increased. The annual 
budget for repairables takes approximately 60-70 percent of the total annual stock fund 
budget of the Korean Air Force, however the essentiality of high quality repairables in- 
ventor)’ management was not realized. [Ref. 1: p. 1] 

In this chapter, the Korean Air Force Logistics system will be described. The 
chapter specifically deals with repairables management, organization, system parame- 
ters, and mathematical models related to repairable item management. This chapter will 
provide the basis for the analysis and suggested improvements to the Korean Air Force 
repairables inventor)-’ management system which will be described later. 

B. THE KOREAN AIR FORCE LOGISTICS COMMAND ORGANIZATION 

The Korean Air Force Logistics Command (KAFLC) is responsible for managing 
the allocation of the logistics resources to support maximum weapon system operational 
availability. It functions as an intermediate echelon command among the tactical units 
of the Korean Air Force under the policies and planning of Air Force Headquarters. 
Figure 1 shows the organization of the KAFLC. The following paragraphs will explain 
the basic missions and responsibilities of the six major divisions of the command. 

The DMM (Directorate of Materials Management) is the focal point of material 
management for the Korean Air Force. It procures all materials according to its esti- 
mates of requirements, distributes the material to all of the tactical units and the sup- 
porting units under the policies of the KAFLC. Thus, the DMM is the equivalent of a 
wholesale level Inventory Control Point (ICP). The DMM manages approximately 
210,000 items which are allocated among the hundreds of item managers. To manage 
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Tigurc 1. The Organization of tlic Korean Air Force Logistics Connnand 

cflcciivcly, DMM relies heavily on the Electronic Data Processing Center's (1;DPC) 
computer services to access historical transaction data, compute demand forecasts, cur- 
rent asset data. etc. 

The DST (Depot of Storage and Transportation) is the centralized warehouse for 
the Korean Air Force where all procured and repaired materials, including consumable'-, 
arc stored, fhey control the physical inventory of the items needed by the ultimate user, 
'fhe DST also deals with the disposal of salvage, scrap, c.xccss and obsolete material. 
Tile other important role of the DS’f is transporting all materials to arrive at the tactical 
unit when required. 

Ihe D.VIE (Depot of iVlaintcnancc and Equipment) is the iii-house depot level 
maintenance organization for the Korean Air Force. The primaiw task of the DME is 
the overhaul of aircraft on a scheduled basis (preventive maintenance). It performs the 
maintenance of aviation end items which arc beyond the capabilities of the the interme- 
diate and organizational maintenance level facilities. These items would be the major 
end items of an aircraft such as the engine, fuselage, gcarbo.v and the direct ground 
support equipment for the operation of the aircraft. 
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The DMEC (Depot of Maintenance, Electronics and Communication) performs the 
maintenance of radar equipment and its subassemblies, ground communication equip- 
ment and weather forecasting equipment. The DMEC also performs preventive main- 
tenance. 

The DMA (Depot of Maintenance and Ammunition) performs the maintenance of 
the precision measurement equipment (PME), airborne equipment and armament of the 
aircraft, such as laser guidance system. 

The EDPC (Electronic Data Processing Center) provides logistics software devel- 
opment and support, establishes job processsing standardization, collection of data re- 
lated to preventive and corrective maintenance, provides analysis and guidance for each 
of the other major divisions. The major objective of the EDPC is to assist each division 
in performing their specialized logistics missions within the Korean Air Force Logistics 
Command. 

C. REPAIRABLE ITEMS MANAGEMENT SYSTEM 

1. The General System Overvie\v 

In the Korean Air Force an item of supply is designated as a repairable if it can 
be repaired faster and less e.xpensively than it can be procured. Weapon systems in- 
stalled in aircraft have become increasingly sophisticated and complex. Hence, many 
weapon systems are made up of a number of subsystems which in turn are comprised 
of several repairable modules. Often the complexity of these individual modules are such 
that the personnel and equipment are not available at the end user level to repair failed 
units. Consequently, these modules are designated as repairables and failed units are 
returned to designated repair activities for repair. Therefore, repairables management 
has become an essential part of the Korean Air Force supply system in terms of decision 
making. Figure 2 shows the current Korean Air Force repairables management cycle. 

The Korean Air Force has two kinds of maintenance operations for repairables, 
the base maintenance (organizational and intermediate level) and the depot level main- 
tenance (DME, DMEC, DMA). It should be noted that not all depot level repairables 
are repaired at the Korean Air Force facilities. Because of the complexity, some repairs 
are done by the USAF under Foreign Military Sales (FMS) using commercial contrac- 
tors. As depicted in Figure 2, three kinds of repairables flows are shown. The first is the 
flow of carcasses to base level maintenance. The second is the flow of carcasses to depot 
level maintenance. The third level is the flow of carcasses from the DST to the USAF 
(the activities that coordinate with USAFLC for repair and procurement follow FMS 
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Figure 2. Tlic Korean Air Force Repairable Management Cycle 
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procedures). Note that the third flow implies the exchange of a carcass for a new item 
from the USAFLC. 

When a failure of a system or subsystem occurs, the equipment is inspected to 
determine the cause of the failure. After isolation of the failure to a component or major 
assembly, further inspection is done to determine if base level repair is possible. In case 
of the failure of a major end item, an RFI (Ready-for-Issue) item from the base supply 
is used to replace the failed item and the aircraft returns to operational status. This in- 
sures the operational availability of an aircraft. If possible, the failed item is repaired 
by the Base Maintenance Squadron (organization or intermediate level of repair). After 
such a repair, the item is released to the Base Supply Squadron in serviceable condition. 

For the determination of the maintenance level, the Korean Air Force has a 
special code assigned to each replaceable component of a major end item. Thus, the 
level of maintenance for each component is predetermined by this code. The code is as 
follows. (Ref 3; p. 7] 

First code position 

X expendable item 

N repairable item 

Second code position 

D depot level maintenance item 

F field level maintenance item (intermediate) 

B base level maintenance item 

Third code position 

1 depot level maintenance 

2 intermediate lev'el maintenance 

3 no maintenance action 

After the Base Supply Squadron has issued a serviceable item from its stocks 
and when base level maintenance has been determined to be impossible, the Base Supply 
Squadron sends the failed item to the DST. Of course, all these actions are reported to 
and directed centrally by the DMM. The failed items from each base are placed in the 
DST facility and await repair. These carcasses are turned over to the depot maintenance 
facility (DME, DM EC, DMA) in batches under the approval of the DM.M. In other 
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words, all carcasses from each base are stocked in the DST facility during the quarter. 
.At the end of each quarter, all carcasses are sent to the depot maintenance facility 
(D.ME. DMEC. D.MA). However, at this time, item managers of the DMM consider 
the capacity limit of the depot maintenance facility.l 

2. Operating Characteristics of Repairables Management 

As mentioned earlier, there are sev’eral depot maintenance organizations in the 
Korean Air Force. However, their concept of operation and relations with the D.MM 
are identical for the DME, D.MEC and DMA. So, only the operation of DME will be 
covered here. The maintenance organization of the Korean Air Force is restricted in 
terms of its capacity to process incoming repairs and in terms of the level of technology 
to deal with the repair of complex systems. Thus, the .MRS (Material Repair Schedule) 
and MRRL (Material Repair Return List) are established to manage these restrictions. 

a. Material Repair Schedule (MRS) 

The MRS is the repair schedule for the depot maintenance facility. It de- 
scribes the type and quantity of repairable material which can be repaired during the 
next fiscal year by repair depots operated by the Korean Air Force. In September of the 
fiscal year, all item managers estimate the demand for each repairable item for the next 
year, and the total repair quantity, based on data from the last several years. The inte- 
grated results for each item are provided to the D.ME, DMEC, and DMA. Then, the 
DME, DMEC and DM.A set up their repair schedules after discussion and coordination 
with D.M.M. The forecasted maintenance requirements which are not covered by the 
MRS are turned over to the MRRL. The repair quantity for the MRS is set up on a 
quarterly basis and carcasses are issued from the DST based on capacity of the D.ME. 

b. Material Repair Return List (MRRL) 

This is the list of repairables for which the Korean Air Force does not have 
depot level maintenance capability, either due to technology or capacity restrictions. In 
the case of an .MRRL repair, the DST sends the carcasses to a USAF depot facility. 
Upon the receipt of the carcass from the Korean Air Force, the L’SAF returns a serv- 
iceable unit to the DST. Serviceable units from both the USAF and Korean Air Force 
depot facility are integrated at the DST to make up the serviceable stocks w'hich are 
available to support for the tactical unit. 

1 This capacity limit established based on the .MRS, see paragraph 2.a. above. 
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c. The Measure of Performance 

The Korean Air Force supply system (wholesale level) uses two measures 
of performance. The first one is fill rate. It is also used by the base level supply system. 
The second measure is supply response time. It measures the length of time clasped for 
base backorders to be satisfied by the stock from the DST. In other words, it is the 
length of time from the placement of an order with the DMM to the receipt of the or- 
dered item. Notice that the DMM uses this more as a priority rule for issuing the stock 
than as a general measure of performance for management of the repairable items. 
Currently, the Korean Air Force assumes that any backorders for spares result in an 
aircraft being "not operationally ready because of supply" (NORS). Thus, the Korean 
Air Force has developed codes w'hich indicate the NORS condition. Each code provides 
the maximum supply response time requirement to fill backorders. Table 1 summarizes 
the supply response time requirements for each demand with priority and Table 2 de- 
scribes the meaning of each NORS code. [Ref 4: p. 81] 



Table 1. SUPPLY RESPONSE TIME REQUIREMENTS 



Priority 


Type of order 


Response time (day) 


Circumstances 


03 


Express 


3 


G, K NORS 


06 


Semi-express 


14 


A, F NORS 


13 


Routine 


30 


Routine 



Table 2. DESCRIPTION OF NORS CODE 



NORS Code 


Description 


G 


Aircraft is in a Totally inoperational condition 


K 


Radar malfunctioning 


F 


Operational, incapable of mission 


A 


NORS is anticipated 
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D. MATHEMATICAL BACKGROUND OF REPAIRABLES MANAGEMENT 

1. DDR (Daily Demand Rate) 

This is the average number of demands per day and it is computed following 
formula:- 



365 

DDR = y^ 



A 

365 



where: 

D, = demand on the ith day of the year 

2. DRP (Depot Repair Percent) 

It is computed in the following manner for each repairable item from the last 
three years of historical data: 



DRP 



RTS 

RTS + .\RTS + COND 



X 100 



where: 

RTS = repair this station 
NRTS = non-repair this station 
CON'D = condemned 

RTS refers to the number of units repaired at the Korean Air Force depot. 
Conversely NRTS refers to the number of units repaired by the commercial contractors 
or by the USAF under FMS policies. Thus, NDRP (Non-Depot Repair Percentage) is 
the complement of DRP. 

3. RCT (Repair Cycle Time) 

The Korean Air Force has two kinds of repair cycle time. The one is the repair 
cycle time associated with base maintenance and the other is associated with the depot 
level maintenance. Usually RCT refers to the time allowed for the depot level mainte- 
nance only. This time allowance is the policy of the Korean Air Force. So, RCT is 
constrained3 to be not less than 30 days and not more than 120 days. 

2 This is computed by the end of year. The time period of this computation starts from the 
fust day of the year to the end of the year. 

3 This is the policy of the Korean Air Force, not a value which is computed from actual 
observations which are recorded in the inventory management database. 



11 



4. RCQ (Repair Cycle Quantity) 

This refers to the quantity of the repairable item which must be stocked to meet 
demands during the repair cycle time. RCQ is applied to MRS items only and computed 
by the following fonnula: 



RCQ = DDR X DRP x RCT 



where: 

DDR = daily demand rate 
DRP = depot repair percentage 
RCT = repair cycle time (days) 

5. OLQ (Operational Level Quantity) 

This is the stockage quantity which is required to cover forecasted demand 
during the time an MRRL item is in the hands of the USAF for repair. OLQ is related 
to MRRL as RCQ is to an MRS item. The formula for OLQ is as follow: 

OLQ = DDR X 60 day 



where: 

DDR = daily demand rate 

6. OST (Order and Shipping Time) 

OST is defined as the time elasped from the initiation of a procurement or re- 
pair order, to the receipt of the order. In the case of a procurement, the Korean Air 
Force constrains the OST to be not less than 120 days and not greater than 365 days. 
For MRRL repairables, OST is constrained to be not less than 220 days and not greater 
than 465 days. In both cases, the upper and lower intervals are adopted to avoid an 
inventory stockout. The 100 days increment in the OST is due to the additional trans- 
portation time for MRRL items from Korea to the continental U.S.A. 

7. OSTQ (Order and Shipping Time Quantity) 

OSTQ refers to the quantity of repairable items needed to meet the demand 
during the order and shipping time. The formula for OSTQ is as follows: 

for MRS items, OSTQ = DDR x NDRP x OST 
for MRRL items, OSTQ = DDR x OST 



uhere: 
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DDR = daily demand rate 

NDRP = non-depot repair percentage 

OST = order and shipping time 

8. SLQ (Safety Level Quantity) 

This refers to the quantity of repairables stocked to meet demand during delays 
in shipping, delays in maintenance, or unexpected increases in demand. The formula for 
SLQ is as follows; 



for MRS items, SLQ = ^/3 x {RCQ -I- OSTQ) 
for MRRL items, SLQ = x OSTQ 

where: 

RCQ = repair cycle quantity 

OSTQ = order and shipping time quantity 

9. RO (Requisition Objective) 

Even with the assistance of the existing computer system, it is hard for the item 
managers to manage the more than 2,000 items (per item manager) that are assigned to 
them. Thus, the DMM designates items that have two or more requisitions per year as 
requisition objective items. The item managers will monitor these items more closely. 

Currently, the DMM applies a periodic inventor}’ review model which is based 
on a policy of reviewing and repairing at fixed regular intervals to bring inventor}’ levels 
back up to the Requisition Objective (RO). Repair orders are placed at predetermined 
interv’als■^ (the expected demand between intervals plus some allowance for the variabil- 
ity of demand must be considered). 

The requisitioning objective consists of three terms. These are the order and 
shipping time quantity (OSTQ), the repair cycle quantity (RCQ) and safety level quan- 
tity (SLQ). The requisitioning objective is the sum of these components. In Figure 3. 
the repair order quantity is simply the difference betw’een RO and the on hand inventor}’ 
at time of review’. 

Since the system parameters such as OST, RCT, DDR, DRP and NDRP differ 
at each review’ period, the level of the RO may also be different at the time of each re- 

4 This predetermined interv’al is the repair induction inter\'al on the MRS. The current length 
of the interval is 2 months. 
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Figure 3. Tlie Korean Air Force Repairablcs Inventory Periodic Revieu 

view. Inventory reviews are usually done every 2 months. The formula for RO is shown 
below: 

for MRS items, RO = SLO + RCO + OSW 

where: 

SLQ = V'3 X (RCQ + OSTQ) 

RCQ = DDR X DRP x RCT 
OSTQ = DDR X NDRP x OST 

for MRRL items, RO = SLQ + OLQ + OSTQ 



where: 

SLQ = 73 X OSTQ 
OLQ = DDR X 60 days 
OSTQ = DDR X OST 
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10. History of RO Computation 

This section shows the history of the RO computation formula. The Korean 
Air Force tries to determine optimal stock level for the repairables inventory manage- 
ment system. However, the results of changing the RO formulation without adequate 
consideration of the system parameters may result in excess stock for some items or 
stockouts for some items. 

a. Before 1973 



RO = OLQ + SLQ + OSTQ 



where: 

OLQ = DDR X 75 
SLQ = DDR X 15 
OSTQ = DDR X 30 

b. Between 1973 and 1976 



For MRS, RO = DDR x 120 

For MRRL, RO = RCQ + SLQ + OSTQ 



where: 

RCQ = DDR X 15 

SLQ = ^ 3 X (RCQ + OSTQ) 

OSTQ = DDR X 84 

c. Between 1976 and 1979 

RO = RCQ 4- SLQ 4- A’CD 4- OSTQ 



where: 

RCQ = DDR X DRP x RCT 
SLQ = 3(RCQ 4 - OSTQ 4 - NCQ) 
NCQ = DDR X NDRP x NCT 
OSTQ = DDR X NDRP x OST 
NCT = NRTS Condemn Time5 



5 This (NCT) is a time spent for decisioning a condemnation among NRTS items. 
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d. Between 1979 and 1984 



RO = RCQ + SLQ + OSTO 



where: 

RCQ = DDR X DRP x RCT 

SLQ = (DDR X NDRP x POST) + RCQ 

OSTQ = DDR X NDRP x ROST 

e. After 1984 

For MRS, RO = RCQ + SLQ + OSTQ 

where: 

RCQ = DDR X DRP x RCT 
SLQ = ^3 X (RCQ + OSTQ) 

OSTQ = DDR X NDRP x OST 

For MRRL, RO = OLQ + SLQ + OSTQ 



where: 

OLQ = DDR X 60 
SLQ = V'3 OSTQ 
OSTQ = DDR X OST 

E. SUMMARY AND PROBLEM STATEMENT 

As mentioned earlier, militar>’ inventory managers tend to want to increase their 
operational availability rather than minimize inventory costs. However, in the case of 
the Korean Air Force repairables management process, the mathematical model is in- 
suflicient to adequately model the stochastic characteristics of the inventory manage- 
ment problem. In the RO (Requisition Objective) computation it does not specify the 
wear-out rate or the regeneration rate, which are important parameters in dealing with 
repairables inventory management. 

The DDR (Daily Demand Rate) computation is the mean of a stochastic process. 
However, the current inventor^’ model treats demand as if it were deterministic. In the 
formula for SLQ (Safety Stock. Quantity), they considered the stochastic nature of de- 
mand but the actual probability distribution is not considered at all. As evidence of the 
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model's shortcomings, the Korean Air Force repairables management tends to stock 
large quantities of spare parts and components in order to avoid stockout situations. 
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111. OVERVIEW OF REPAIRABLE MANAGEMENT MODEL 



A. INTRODUCTION 

Inventory systems such as the one described in the previous chapter can be classified 
rather broadly as a multiechelon system with repair. Analytical studies of multiechelon 
systems have shown the computation in such models to be hard, so either simplifying 
assumptions or approximations are necessary. The introduction of the repair aspect into 
the system certainly complicates the structure of a model even at single echelon levels. 

Suppose a repairable system, consisting of one Inventor}' Control Point (ICP), one 
stock point and one overhaul and repair activity controls the inventory' of a single item. 
Demands from various tactical units are placed only at the stock point. The system has 
continuous updating of records, i.e., transaction reporting. When items wear out or fail, 
the tactical unit can either scrap the item or return it to the depot repair facility (D.ME, 
D.M.A, D.MEC). After inspection, the depot can either scrap the item or repair and re- 
turn it to the RFl inventory. 

In this chapter, we will discuss various repairable inventory' models based on two 
distinguishing characteristics. One is a deterministic model, the other is a stochastic 
model. In case of deterministic model, we will describe two models, "continuous sup- 
plement" model and "substitution" model. For the stochastic case, we will describe five 
models based on current usage in U.S military'. Each description will include model 
formulation and solution techniques. 

B. DETERMINISTIC INVENTORY MODEL 

These deterministic inventory models were developed by David A. Schrady and 
others between 1966 and 1967. These deterministic inventory models are described in a 
research report for the U.S Navy repairable inventory management system [Ref 5, 6, 
7]. The first model calls for the repair activity to induct a batch of carcasses when the 
NRFl inventory level reaches a certain point. With a deterministic system, this rule in- 
sures regularly spaced repair inductions of fixed batch size, which is a simple process for 
the repair activity to manage. Repair batch sizes can be maintained if there is contin- 
uous supplementing of repaired RFI items with new procurement. The procurement 
order will be issued with enough lead time so that when the on-hand RFl stock drops 
to zero, a procurement quantity will arrive. The two inventories, RFl and NRFI, will 
have time histories as shown in Figure 4. Note that in this, the "continuous supplement" 
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model, the procurement trigger is in the RFI inventory and the repair trigger is in the 
NRFI inventor}'. [Ref. 6, 7; pp. 9-10, pp. 392-393] 

The second model is suggested by noticing that there is a trade-olT between stock 
held in RFI condition and stock held in NRFI condition. NRFI carcasses are a legiti- 
mate resource because they require only repair to be restored to full usefulness. But the 
value of this resource is less than the value of the RFI resource by at least the cost of 
labor and replacement parts. Thus, if inventor}' is to be held in the system, it would be 
slightly cheaper to hold it in NRFI condition than in RFI condition. [Ref 6, 7: p. 10, 
p. 393] 

Notice that the "substitution" policy minimizes the RFI inventory while the "con- 
tinuous supplement" policy minimizes the NRFI inventor}' [Ref 6: p. 13]. The ne.xt 
section describes the mathematics of these two models. 

1. Analysis of The Continuous Supplement Policy Model 
a. Notation 

Under the assumption of deterministic demand and lead times, it is not 
necessary to maintain a safety stock. Whenever the RFI inventory reaches zero, a pro- 
curement quantity will arrive. Thus, the system is never out of stock. Between pro- 
curement arrivals, depleting RFI stock is replenished with repaired items. Also, a 
repetitive system will be established regardless of the initial provisioning policy, so that 
it is sufficient to analyze only one cycle to determine system characteristics. Figure 4, 
shows basic characteristics of the "continuous supplement" policy model. The basic 
notation as follows [Ref 5: pp. 12-22]: 

D = annual demand rate 

r, = procurement cycle 

r, = repair cycle 

t, = procurement lead time 

tj = repair lead time 

To = field recover}' rate 

/•j = overhaul and repair recover}' rate 

Q, = procured lot size 

Qi = repaired lot size 

h^ = holding cost for RFI 

hi = holding cost for NRFI 

A\ = reorder point for procurement 
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X2 — reorder point for repair 

n = number of repair cycles per procurement cycle, n = Til 
Ax = fixed procurement cost per one cycle 
A2 = fixed repair cost per induction 

b. Model Formulation 

The objective of the model is to determine an optimal procurement quan- 
tity Q\ and an optimal repair quantity that minimize the average annual cost. These 
optimal values, coupled with respective reorder points, A'l and X^, constitute an optimal 
operating doctrine. To determine the average annual variable cost, the costs per repair 
reorder cycle must be investigated. The product of per cycle costs and the number of 
cycles per year will then yield the average annual cost. The total cost for any given cycle 
Tx is the sum of procurement, repair and holding costs. 

Since there is only one procurement per reorder cycle, the variable order 
cost is the actual cost of the items ordered and can be expressed as QxCx- The fixed 
procurement cost for one cycle is Ax as mentioned. To determine the repair costs per 
cycle it is necessary to compute the number of repair cycles per procurement cycle ( 
n = Tx I Tj, n is an integer). In Figure 4 : 



Tj = r, -b 



where; 

tx = KT2, 0<K<1 

h = (l - Aor, 

The cost of items repaired per procurement cycle is C2Q2 and the fixed repair cost is 
/Ijrt, where A2 is defined to be the fixed repair cost per induction. The holding costs per 
procurement cycle will be hxAj^ -h /i2/ir2' where Ar^ is the area under RFI curve and 
Aj^ is the area under the NRFl curve. To compute the area under RFI curve, consider 
Figure 5 showing RFI inventor}' for one procurement cycle. Since D is constant and 
known, the area of U, A^, is given by: 

Au ~ ~2 



20 







SCK MP 



mi 





T. 



-H 



p.rt 






1 




Figure 4. Repair Cycle and System for "Continuous Supplement" Model 
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where: 



h = kt, 

a = Q^ — D ti 
The formula will reduce to: 



Au= (h2i - Di]) I 2 

Since n is an integer, the area of V is the sum ofn-1 trapezoids each having 
one side of length Tj. To determine the area of each trapezoid, a recursive relation is 
developed. Let: 

Ay. = area of the ith trapezoid in V (see Figure 4) 

= Ay^ + Ay^ + + Ay^_^ 

where: 

Ay, =Y U2Q, + 2{Q, - t,D)~ DT,-\ 

Ay, = Y T,\_AQ, + 2{0,-t,D)- 3D7J 
therefore, the general formula for the area of V: 

Ay^ = ^ U2iQ, + 2{Q, - h D) - (2/- l)DT,-} i= 1,2,....,a7 - I. 

n-l 

So the area of V is, Ay= "^Ay^ 



Ay= -T U2Q,Zi + 2{n-\)iQ, - t,D) - DT,Y.{2i-\y] 

^ /= 1 i= 1 

= Y T,iQ,n{n-\) + 2{n-\)Q, - 2{n - \)t,D - DTA^ ~ \) + DUn - \)] 
= ^ r,[2(«- l)(a - hD) + {n^ - n)Q, - («- 1)^Z)7',] 

Finally, area of W in Figure 4, A^ is easily calculated as follow; 
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where: 



r, = (1 - k)T, 

b = nO. + (2i - hD) - (« - \)DT^ 

Thus, the area of W; 

Aw = jhinQ2 + Q^ - hD - {n- \)DT{\ 
The total area under RFI curve is: 



Aj — Au + Ay + Ayr 



= r,0, - \di\ + ^ T,{2{n - 1)(0,) - Di, + {n^ - n)Q, - {n^ - n)DT,] 

Since n is an integer and the build up rate of NRFl items, r^D, is constant, 
the area under the NRFl curve is simply n times the area under the repair cycle curve 
is TjO, / 2/-;. So the NRFl holding cost per procurement cycle is / Ir-^. The total 

cost per procurement cycle becomes: 



^r, = Q\C\ + A^ + 



Ti 



C2Q2T, 

T2 



+ 



lhQ2h 



+ hiA-p^ 



The total averaae annual cost is then: 



K = 



T, 



giC, 

T, 



+ 



T, 



+ 



T, 



■ + 



C2Q2 



+ 



^2Qi 



+ 



~Tr 



Simplifying the total average annual cost: 



K = 



A^RD 



^202 ^\Q\ 

C,RD + A,,,,„D + + -if- 



h\Q2 

— h^kQ2 + 2 



where: 

7*2 = Q 2 I V 2 D 

TJ T 2 = n , n = 1,2,3,.... 

r, = A: Tj = kQ2 / 0 < A < 1 

r, + = 7*2 

Ti = QJ \ - ^20 = QJ RD 



23 



y? = 1 - >'0^2 

Note that the terms C^RD and are independent of and Q^, and hence do not 

affect the operating doctrine. Therefore it is appropriate to redefine the average annual 
cost of ordering, repairing and holding as follow: 



K = 



A^RD 



AjYQrjD 

Q2 



2'2 



+ 



Mi 

2 



h\kQ2 + 



^1^2 

2 



From differential calculus the optimal Q^ and will be those that satisfy the equations, 



cK 

oQ^ 



= 0 , 



dK 

0Q2 



= 0 



Taking the partial derivatives results in: 



dK _ A-^RD ^ 

ce. ■ Ql ^ 



cK 

^Qi 



Ql 2 ,, 




Solving the equation to get the optimal procurement and repair quantities, 



. ^ l 2A,RD . ^ I 2A,vlD 

V /'I ’ V *2 + '■ 2 * 1(1 - 2*) 

2. The Substitution Policy Model 
a. Notation 

Notice that this model minimizes the RFI inventory. Figure 5 shows basic 
characteristics of the "substitution " policy model. The basic notation as follows [Ref 
6, 7: pp. 13-17, pp. 393-396): 

Qp = procurement quantity 
Qp = repair batch size 
d = demand rate 

(1 — r) = scrap rate, r is the recover}' rate measured 
ip, tp = procurement and repair lead times 
Ap = fixed procurement cost per order 
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.-Jr = fixed repair batch induction cost per batch 
/:, = RFI holding cost 
/i, = NRFI holding cost 
T = system cycle time 

b. Model Formulation 

Referring to Figure 5, define 7, as the time period during which inductions 
are suspended and the repair is simply accumulating NRFI carcasses. It may be deter- 
mined that: 



7 ; = {Qp + QpMd 

Next, let n be the number of inductions per cycle. In general, it will not be possible to 
insure that the last induction before begins will reduce the NRFI stock to zero. The 
residual NRFI stock when T, begins will be some fraction of the net loss in NRFI items 
per repair cycle, where the repair cycle is the time between regular consecutive in- 
ductions. The net loss of NRFI inventoiy over the period between successive in- 
ductions; 



^d{Qp I d) — Q[i — — Qr{\ — r) 



Thus, the residual when T„ begins will be some fraction of Q^{\ — r), call it — r) 

where 0 < JS < 1. While the /? factor is unavoidable due to the requirement that there be 
an integral number of repair batches in each cycle, we shall set ^ = 0 in subsequent de- 
velopments. 

Now the first batch after inductions are resumed takes the amount from 
the NRFI inventory. Subsequent inductions cause a net reduction of only Q^(l — r) 
items. Thus, the amount of NRFI items available for the (n-1) inductions are 
rdT^ — or Qn{r — 1) -f rOp. Dividing the last expression by Qp{\ — r) yield: 



(«-l) 



Qr{^ ~ 0 + ^Qp 
Qpi^-r) 



- 1 + {j^)QpQr 



or 



n = ( 



r Qp 
Qr 
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Figure 5. Repair Cycles ami Sy.slein for "Substitution" Model 
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The system cycle time is given by: 

_ riQR + Qp _ Qp 

~ d ~ {\-r)d 

where: 

^ ~ ( j _ ^ )QfQr- 

From the Figure 7, the area under the RFl curve, over one complete cycle Ai, is then: 
Referring again to Figure 7 the area the NRFI curve: 

‘^2 = + ej] 



Now we get total cost per unit times expression: 



TC _ , d[ 1 - r) Aprd 
T ^ Qp Qr 



hr 



^[.Qr + ( 



)Qpl + hrLQp + Cr] 



where: 



T = 



Qp + riOp 

d 



Qp 

di\-r) 



The total cost under substitution model: 

TC = Ap + fiAr + j ^ ^ {.QpQr + r ^ "i 1 p ^QpQr Qp^ 

The optimal order quantities are obtained by setting the partial deriatives, 
with respect to Qp and Qp. 



8TC ^ ^ ^ Q 



6Qp 



^Qr 



Next we get new formula for optimal order quantities: 
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dTC ^ ->') ^(1 ~r) _v 

8Qp Ql ^ ^ 

cTC ^ _ ^R''d fhr hjY 

8Qr Ql 2 2 

The optimal order quantities for repair and procurement: 



Qp — 



2Apd{\ -r) 
A,(l - r) + V 



= 



2A^d 

/ i , + /?2 



3. Summary and Evaluation 

An important result of deterministic inventory model analysis is the realization 
that economic order quantities do not vary directly with demand but as the square root 
of demand. This nonlinear relationship could explain why many organizations experi- 
ence inventory problems. Intuitive inventory policies are undermined by difficult-to- 
grasp nonlinearities. [Ref 8: p. 141] 

Although the models described in this section are evidently applicable to re- 
pairables management, they are fraught with many limitations. The reasonableness of 
model assumptions and their sensitivity to actual conditions determines the utility of any 
particular model. For example, the two models assume that the demand for a repairable 
item is known with certainty. Also stockouts either are not permitted to exist or are 
backordered and satisfied when replenishments are received. 

In terms of the Korean Air Force repairable management system, the models 
described here are difficult to adapt to the Korean Air Force system. This is because 
of demand variability, relatively long procurement lead time, repair capacity and capa- 
bility, a number of backorder situations, and current system management concept. The 
other important point is that the Korean Air Force often encounters NORS conditions 
which results in decreased operational availability. However, these two deterministic 
models do not deal with stockouts. These reasons didn't allow to use for the Korean 
Air Force repairable management system. 

C. STOCHASTIC INVENTORY MODELS 

As the weapon systems installed in modern aircraft become increasingly sophisti- 
cated and complex, repairable items represents an important subset of the total inven- 
tory of items which are managed by the military supply system. 
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In general, repairable items are supported by a two-echelon inventor}’ and repair 
system. When a repairable item fails at the base level, it is returned to base supply and 
a new serviceable unit is issued from e.xtra RFI stock of the base. If possible, the failed 
item is then repaired by the base maintenance organization and then stored in RFI 
conditions at base supply. Sometimes, however, the failed item must be returned to the 
repair depot where more sophisticated equipment and specialized skills are available to 
repair it. 

When we consider the condemnation of repairables, there should be an inflow of 
newly manufactured items to depot level supply should look at the material flow of all 
echelons, and then, place a replenishment order to resupply the condemned repairables. 
.A.t times, forward base locations will be supplied from another closely located base. This 
is called lateral resupply. For other items, a manufacturer may provide both the source 
of procurement for new assets and the source of repair for failed items. 

The decisions concerning which maintenance levels will repair the failed items sub- 
sequently affects the supply support provided at the organizational level. If an item is 
repaired at the organizational level, repair equipment and maintenance personnel must 
be made available at that level. However, if the item is not repairable at that level, then 
the question is whether the item can be replaced at the organizational level. Organiza- 
tional level repairables will normally be transferred to the next higher echelon of repair 
if the repairs cannot be accomplished at the organizational level. Similarly, repairables 
which cannot be repaired by the intermediate level are usually sent to the depot level for 
repair. 

The difference between stochastic models and deterministic models is that if de- 
mands or future requirements are uncertain, then regardless of the stockage policy 
adopted there is generally some probability that available stock levels will be insufficient 
to meet demand (Ref 9; pp. 187-188]. In this section, we will describe five stochastic 
models based on current usage. 

1. Model Analysis 

a. Analysis of METRIC Model 

Background. The METRIC model was developed over a period of 
years by a research group at the RAND corporation and presented in the literature by 
Sherbrooke (1968) and extened by Muckstadt (1973). METRIC was developed with the 
ultimate goal of implementation and a slightly modified version of METRIC was actu- 
ally implemented by the L’SAF. (Ref 10. II, 12] 
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METRIC is a mathematical model of a base-depot supply system in 
which item demand is assumed to be compound Poisson with a mean value estimated 
by a Bayesian procedure. The model is also capable of determining base and depot stock 
levels for stock levels of a group of recoverable items. Its governing purpose is to opti- 
mize system performance for specified levels of system investment. It is designed for 
application at the weapon-system level, where a particular line item may be demanded 
at several bases and the bases are supported by one central depot. The support depot 
may vary by item as in the item-manager system or it may be fixed as in the weapon- 
system storage site concept. The general purposes of this model are optimization, re- 
distribution and evaluation. [Ref. 13; p. 2] 

The basic METRIC model considers a two-echelon system in which 
independent bases (lower echelon) are supported by a repair depot (upper echelon). At 
the occurrence of a failure (it could be more than one), the failed item is either replaced 
by available base stock or back ordered if the base stock is not available. The item is 
inspected to determine the extent of repair required. If the repair can be made at the 
base, the unrepaired item enters base repair. If the item cannot be repaired at the base 
level, it is shipped to the depot. Upon shipping the item to the depot, the base places 
an order with the depot for a replacement, so that the inventory' position for item i at 
base j can be maintained. 

(2) M athematical A ssumptions 

1. System Objective of Minimizing the Expected Number of Backorders: The objective 
will be to minimize the sum of expected backorders on all recoverable items at all 
bases pertinent to a specific weapon system. Thus, unless all bases are identical, 
the expected number of backorders will vary by base. Expected backorders take a 
fixed period of time and add together the number of days on which any unit of any 
item at any base is backordered. In order to minimize the expected number of 
backorders one divides this number by the length of the period and taking the ex- 
pected value. This yields a number that is independent of period length. This is 
the value we seek to minimize. [Ref. 14: pp. 126-131] 

2. Compound Poisson Demand: We assume that demand for each item is described by 
a logarithmic Poisson process. Compound Poisson processes are generalizations 
of Poisson processes; the compound processes allow the flexibility of incorporating 
more parameters, yet retain the simple analytical properties of the Poisson. The 
logarithmtic Poisson is obtained by considering batches of demands where the 
number of batches follow's a Poisson process and the number of demands per batch 
has a logarithmic distribution. 

3. Demand is Stationary over the Prediction Period: It is assumed that the distribution 
of demand over some future period of interest is stationary'. 
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4. Repair Decision Deponds on the Complexity: The assumption is that the decision to 
repair a unit at base level or send it to the depot is a function only of the type of 
mallunction and the base maintenance capability. 

5. Lateral Resupply is Ignored: When a unit is shipped from base to depot for repair, 
a serviceable replacement will be resupplied from the depot if possible. If the depot 
has no unit on the shelf, the base must wait until a unit emerges from depot repair. 

6. System is Conservative: Consider a particular stock item demanded at base j. We 
assume that a unit of stock has a probability of being repairable at base j, with 
a repair time drawn at random from a base repair distribution with mean A', a 
probability (1 — r^) of being depot repairable, with an order and shipping-time 
distribution having mean 0, and a depot repair distribution having mean D. This 
implies that there are no condemnations or that the system is conservative, as the 
name "recoverable item" suggests. A higher condemnation rate usually indicates 
that the item should be redesigned. The condemnation rate must be considered for 
procurement purposes, but the procurement process is not considered in the 
METRIC optimization. 

7. The Depot Does Not Batch Units of Recoverable Item for Repair: The model as- 
sumes that depot repair begins when the repairable base turn-in arrives at the de- 
pot. This appears to be a reasonable approximation to current depot scheduling 
practice. Since METRIC economizes by buying fewer high-cost items, these are 
the items which are most likely to be in short supply. 

S. Demand Data from Different Bases Can Be Pooled: We assume that demand from 
the several bases can be pooled in some manner so that a composite initial estimate 
of demand per flying hour can be obtained. The pooled estimate can be obtained 
by a simple averaging technique or a more sophisticated procedure such as expo- 
nential smoothing. METRIC multiplies this number by the flying hours per month 
for each item. 

(3; Model Formulation. By mathematical assumptions, demand in this 
system follows a compound Poisson process. A compound Poisson process may be 
thought of as a series of customers who arrive following a Poisson process, each of 
whom can demand an amount that is independently distributed according to a com- 
pounding distribution. Assume that item i is stocked at each of base j, and the cus- 
tomers who place demand for the item at each base have a known mean arrival rate of 

j= 1,2,3, ,J. When a customer arrives at base to place one or several demands, he 

turns in an equal number of carcasses. These carcasses can be repaired at base level with 
probability of r,^, while (1 - is the probability that they must be repaired at the de- 
pot. The arrival of carcasses from base j at the depot is described by a Poisson process 
whose mean is (1 — r,^ times the mean of the Poisson customer arrival process at the 
base j. Therefore, the total demand at the depot for item i is compound Poisson, with 
mean customer arrival rate; 
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y=i 

where: 

/y = mean arrival rate fot item i at base j 

Let fij be the mean demand per customer at base j. Then, the mean 
depot demand rate for unit i is: 

J J 

y=i y=i 



where: 

6^j = the mean rate for item i at base j 

In the special case of the logarithmic Poisson process, the probability 
that X customer demands are in the repair/ supply process is negative-binomial with pa- 
rameters q and K (Note: K = aT I In ^ where ). is mean customer arrival rate and T is 
average resupply time). 



PixlAyT^ = (K+ X - l)\ 



(q - 1 )' 

{K - l)!.v!/+^' 



X = 0,1,2, , q> I, K>0 



where: 

q = the variance to mean ratio of P{x \ AijT/j) 

K= W, I {q - 1) 

Tij = average resupply time for a demand for item i at base j. 
Atj = base repair cycle time 
D = depot repair cycle time 
= the order and shipping time 
<5 (5(3) = average depot delay 

Then: 



Tij = rijAtj -f- (1 — rij){Oij + ^(5<;>)} 
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Since it takes an average of D time units for an arrival to complete 
the repair process, the probability distribution of the number of units in the depot repair 
cycle is compound Poisson with mean of >.Z). Hence, the expected number of units back 
ordered at the depot is: 



50(Soi/£>) = - So)P(.xi;.z)) 

As mentioned in the beginning, the objective of METRIC is to min- 
imize the sum of backorders for all item i and for all bases j within a budget constraint. 
Thus, the .VIETRIC problem will be represented as follows; 

/ j 

minimize 

i=\ j=\ 

I J 

subject to y^CjSij < C 

i=\ j=\ 

Sjj >0, !</</, 0<j<J 



where: 

S,j = the decision variables, stock for item i at base j 

C = the total amount of budget available 

C, = the cost of item i 

S,o = the depot stock for item i 

''4^ Solution Technique. The .METRIC problem is solved by using the 
generalized Lagrange Multiplier method suggested by Fox and Landi [Ref. 15: pp. 
258-261]. Let cP be a Lagrange Multiplier associated with the budget constraint. The 
Lagrange function is written: 

i J 1 J 

/=1 ;=1 i=l /=! 

The auxiliary' problem attempts to minimize this equation. By trial 
and error, we try to find the value of O which satisfies a given constraint. Therefore, 
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we need to solve the above equation for several values of <t>, and choose that value of 
O for which the required resources are closest to the budget limit. Fox and Landi sug- 
gest a binary search procedure. Their computational experience was that at most six 
bisections were required to obtain budget allocations that were within one half of 1% 
of the original budget C. 

The objective function is separable in the items. Dropping the sub- 
script i in the original problem allows us to rewrite the equation as: 

1 

min Y,BOiSj\XjTj) - OQSy - cPQS,o 

/=i 

Since BOj{Si \ AjTj) is discretely convex for a given the optimum 
base level is obtained by simply finding the smallest non-negative integer satisfying: 

BOjiSi + \\?.jTj) - BO{Sj\).jTj)>OqS,o 

(5) Model Evaluation. The current Korean Air Force Logistics system's 
objective is maximizing operational availability subject to a budget constraint. This is 
similar to the METRIC model. The objective of minimizing the expected number of 
backorders is related to the objectives of the Korean Air Force Logistics system. How- 
ever, two assumptions are distinctly violated. One is batch size repair policy is not used, 
the other is lateral resupply is ignored. Currently, these assumptions are not used in the 
Korean Air Force because of relatively long procurement leadtimes and repair capability 
limitations. It may be possible for the Korean Air Force computer resources to handle 
the data requirements for the METRIC model, but it is doubtful for actual usage. 
However, the .METRIC model is obviously related to the Korean Air Force in terms of 
its basic concept and system objective. 
b. Mod-METRIC Model 

( 1) Background. Mod-METRIC model was developed by Muckstadt to 
deal with problems which the .METRIC approach did not consider. Mod-.METRIC 
considers the relations in parts hierarchy and tries to solve this multi-indenture level 
problem. .Mod-.METRIC was implemented by the U.S. Air Force as the method for 
computing repairable stock levels for the F-15 weapon system. [Ref 16: p. 472] 

Most repairables contain subassemblies or other components which 
are also repairable. For example an aircraft engine, it may have modules for intake, 
combustion and exhaust. If an engine fails, it is replaced by a serviceable engine from 
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base stock. The failed component is then repaired either at the base or depot level de- 
pending upon the complexity of repair. A serviceable module from the base stock, if 
available, will replace the failed module, and the repaired engine is placed in base engine 
stock. 

Q) Marhematical Assumptions. The assumptions stated in METRIC are 
also applicable to Mod-METRIC except that Mod-METRIC assumes that the demand 
process is the simple Poisson process. [Ref. 16 pp. 475-476] 

<’3j Model Formulation. In METRIC, the objective is to minimize ex- 
pected backorders for all items subject to an investment constraint. The Mod-METRIC 
objective, however, is to minimize the backorders for the end item, subject to an invest- 
ment constraint on the total dollars allocated to the end item and its components [Ref. 
16, 17; p. 475, p. 39]. This difference is caused by the hierarchical maintenance re- 
lationship between the module and the end item. As an example, an engine backorder 
indicates that an aircraft is missing an engine and is unavailable to perform its flying 
mission. .Modules, on the other hand, are used only to repair engines. A backorder for 
a module only delays the repair of an engine. The impact of module backorders and 
engine backorders is clearly not the same. 

Let Tj denote the average engine resupply time at a base. By the na- 
ture of repairables, depends on several factors. When it is repaired at base level, 7]- 
would be the time it takes to flow through the base maintenance system. If it is repaired 
at the depot, 7, consists of the time to place the depot order for a serviceable part and 
to receive the part from the depot, assuming a sen'iceable asset is on hand at the depot. 
However, when there are no serviceable assets on hand at the depot, an additional delay 
is included in the resupply time: 



7, = + (1 - ri){Ai + 6{S,D)] 



where: 

r, = the probability an engine will be repaired at base i 

Bi = the average resupply time, given an engine is repaired at base i 

A, = the average order and ship time for an engine at base i 

S, = the depot engine stock 

D = the average depot repair time 

SS^D = the delay days per demand 
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SS^D can be derived in the following manner. The expected number 
of engines back ordered at the depot is; 



BO{S,\;.D) = XI (X - S,)P{AV.D) 
x>s. 



where: 

/ = Z(1 - r,)A, 

/, = the daily engine removal rate at base i 

In other words, this expression is the expected number of units on 
which delay is being incurred at a random point in time. Dividing this expression by the 
expected number of demands per day yields a statistic which has the dimension of delay 
days per demand: 

5S,D = B{S, \;.D) I ;. 

= expected backorders, 'expected daily demand 

Further, B, can be divided into two components. That is, 5, is equal 
to the average remove and replace time, given that the necessary^ module is available, 
plus the expected delay due to the unavailability of the module which is required to re- 
pair the engine. Then: 



Bi = Ri -f A,- 



where: 

R, = the average repair time at base i if modules are available 

A,= the average delay in base engine repair due to unavailability of a needed module 

Further, assume that the engine should be repaired by the failure of 
module j. Then, A,y is the expected delay in engine base repair time due to a back order 
on module j at base i. Thus: 






where: 
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= average number of daily removals of module j at base i 

S. = stock level of module j at base i 

T, j = average resupply time for module j at base i 

Then: 



Tij - rijB:j + (1 — rij){Aij + djDj) 



where; 

r,. = the probability that a failure isolated to module j will be repaired at base level 
B,j = average base repair time for module j at base i 

A, = average order and ship time for module j at base i 
D = average depot repair time for module j 

<)= I(X- S,,)/>(A'|0,Z))/0,Z), 

= the stock level of module j at the depot 

B, = h-,{\ - r,) 

Consequently, the expected delay in engine repair at base i due to 
module unavailability is: 



J 

;=i 

Thus, expression 7i represents all the system components including 
the depot engine and modules stock level, and the base module stock level. By using this 
relationship, the engine and module stock levels can be derived. Finally, the math- 
ematical statement of the problem is: 

/ 

min - WA', I 

, = 1 A‘> 5, 

/ J J 

subject to 'y' (QS; + 

/=1 y=I y=l 



^CySy) -I- '^CjS^J -b C[^So<C 



where; 

S, = stock level of base i 
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Q = unit cost of an engine (end item) 

Cj = unit cost of module j 
C = dollar budget limit 

(4) Solution Technique. Unfortunately, the equation above is not sepa- 
rable, because T,j is a complex function of the S,y. Thus, Muckstadt recommended that 
it should be broken down into two parts; the component subproblem and the end item 
subproblem [Ref 16: p. 477]. Even after the break-down of the problems, however, it 
requires the solution of many subproblems each of which corresponds to a particular 
division of the available budget between components and end items. The solution pro- 
cedure is as follows. 

First allocate the budget C into C, and Cj for components and end 
items repectively. Allocate C, among components so as to minimize the expected end 
items repair delays summed over all bases subject to the budget constraint C, . This 
problem is mathematically represented as: 




j I 

subject to ^ C, 

;=i 7=1 

This problem can be solved using the METRIC technique. Given the 
result from this step, compute the average resupply time Tj for the components of each 
base. Then, allocate the remaining budget so as to minimize the expected end item 
base backorders. The METRIC budget allocation procedure may again be used. The 
steps provide a set of proposed stock levels upon a given allocation of the budget among 
the end items and components. These steps then are repeated several times using new 
values for C, and to establish the best allocation. 

(5) Model Evaluation. Basically, the Mod-METRIC model has a similar 
structure to METRIC. All assumptions used in Mod-METRIC model are applicable to 
.METRIC except for the compound Poisson demand. In other words, this model is the 
same as .METRIC for the Korean Air Force repairable management system. The only 
disadvantage is that the .VIod-.METRlC model was developed for new weapon systems 
such as the F-15 aircraft. However, the Korean Air Force has a lot of non-standard 
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weapon systems. These non-standard weapon systems are causing problems such as 
long procurement leadtimes. 

c. Dyna-METRIC Model 

n J Background. Dyna-METRIC was developed by the R.A.ND Corpo- 
ration to provide an analytic method for studying the transient behaviour of 
component-repair and inventory systems under time-dependent operational demands 
and logistics decisions like those that might be experienced in wartime [Ref 14, 18]. 
Note that the past work regarding the repairable item stockage prior to Dyna-METRIC 
only dealt with a steady state inventoiy system with constant average demand and ser- 
vice rates. These steady state assumptions may provide a good approximation during 
peace time operations. But in wartime demands for components may suddenly increase 
ver>' much relative to the previous peacetime operation and then may decrease gradually 
or, in some case, drastically due to attrition in the system. [Ref 14: pp. 1-6] 

A key characteristic of Dyna-METRIC is its ability to deal with the 
dynamic or transient demands placed on component repair and inventory support 
caused by time variables in a scenario that includes sortie rates, mission changes, phased 
arrival of component repair resources, interruptions of transportation, and the like, all 
of which would be experienced in wartime. It computes how given resource levels and 
process times would contribute to the wartime capability. By exploiting the mathemat- 
ical structures of its underlying equations, Dyna-METRIC suggests the alternative cost 
effective repair or stockage resource purchases which would achieve a target aircraft 
availability goal throughout the wartime scenario. 

Dyna-.METRIC considers a three echelon inventoiy repair system. 
Each base has an in-house repair facility which may have various test and repair capa- 
bilities. This base repair facility may be supported by several Centralized Intermediate 
Repair Facilities (ClRFs). Each operating base is capable of doing only limited types 
of maintenance, usually limited to simple removal and replacement operations at the 
flight line. It should be noticed that some of the bases have a CIRF while others do not 
have direct flows of parts to the depot. A depot is represented as existing outside of the 
model. It is seen from the model's point of view, as an infinite source of supply located 
some order and ship time away. 

^2) Mathematical Assumptions. The major assumptions of 
Dyna-.METRIC which distinguish it most from the others are shown below'. 

1. Demand for items are generated by a nonhomogeneous Poisson process with in- 
tensity function m(t), and mean value function r(t)= fm{s) ds. The functions. 
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m(t) and r(t). will be defined later. The t denotes an arbitrary time. Thus, the de- 
mand process is dynamic. 

2. The repair process is independent of the arrival process and has slack repair ca- 
pacity so that each repairable item demanding repair immediately receives services 
with average service time based on the function F(s,t). F(s,t) will be defined later. 

3. No more than one subcomponent can fail or be demanded in the repair of each 
assembly. 

(3) Time Dependent Pipeline Equations. Since the major objective of the 
system is to avoid the loss of aircraft mission capability due to a shortage of correctly 
functioning components on the aircraft, it is necessary to compute the number of com- 
ponents awaiting repair, being repaired, being on the way to and from another echelon 
of repair, and partially repaired but awaiting spare parts. Each state is a pipeline seg- 
ment is which characterized by a delay time that arriving components must spend in the 
pipeline before exiting the segment. The model expands each component's expected 
pipeline size into a complete probability distribution for the number of components 
currently undergoing repair and on order, so the probability distribution for all compo- 
nents can be combined to estimate aircraft availability and sorties. [Ref 14; pp. 7-23] 

Under the assumptions that the probability distribution of repair 
time is independent of the failure process, the average number of components in the re- 
pair pipeline will be: 



dT 



where: 

d = average daily failure rate 
T = average repair time 

With the further assumption that demand has a Poisson probability distribution, the 
probability that there are K components in the pipeline at any point in time would be: 

P{ K in pipeline ) = / K ! 

However, in Dyna-METRIC the demand is a function of time so that: 
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d(t) = (failure per flying hour) 
x(flying hour sorties at time t) 
x(number of sorties day per aircraft at time t) 
x(quantity of the component on the aircraft) 
x(percentage of aircraft with the component) 



Also, in place of a constant average repair time, T, the dynamic model uses the proba- 
bility that a repair started at time s is not completed at time t. That is: 

F(t,s)= Probability(component entering at s is still in repair at t) 

= Probability( repair time > t-s when started at s) 

The average number of components in the pipeline is derived by combining those two 
functions. Consider only those components that arrived in an interval of time, As. cen- 
tered at time s. Then: 



A/.{r.s) = expected number of components in the repair pipeline 
d{s) = daily failure rate at time s 

F{i,s) = Probability that a component is not out of repair by time t 
As = inters’al of time centered at s 

If we assume that the number of failures arriving in the interval As is independent of the 
number of failures arriving in similar intervals centered at other time other than s and 
F(t) is independent of the probability distribution generating the demand rate, then: 



A/(r,s) = d{s) F{i,s)As 



where: 




Further assume that As is very small, so that: 



/(r) = d{s) F{r,s) ds 



0 
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With the additional assumption that the component failure probabil- 
ity distribution is Poisson, /(/) is the mean of a time varying (nonhomogeneous) Poisson 
process. That is, the probability of K components in repair at time t is: 

P{K) = I K ! 



where: 



0 

(4) Time Dependent Component Performance Measure. The component 
measures typically computed by the Dyna-METRIC model are [Ref 14: pp. 24-29]: 

R(t) = ready rate at time t 

FR(t) = fill rate at time t 

EB(t) = expected back orders 

VBD(t) = variance of the backorders 

DT(t)= average cumulative demands by time t 

S(t)= supply level at time t 

The readv rate is siven bv: 

S(l) 

R(t) = Y^P[(K\ >.(:))} 

Since the definition of the fill rate is the probability that a component 
will be available when a demand is placed, it is therefore the probability that demands 
have left at least one component available, that is, the sum of the probabilities of de- 
mands less than the stock level. Expected backorders are given by: 

EB{t)= Y.{K - sit)]P{Kim 

K>S{[) 

S 

= ;.(/) - s(i) + ^(/){5(/) - iQP{Kim 

A'=0 
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For K greater than s(t), there will be backorders of [K — s(r)}. The 
probability of any demand level, that is K, is P{K j /(r)}, and the expected value of the 
backorders is merely the product of the various values the backorders can take on times 
the probability of a demand at that given value. The variance in backorders is given by: 



VB{i) = 

K>S{r) 

(5) Time Dependent Optimal Determination of Spare Parts. The fact 
that pipelines have time-dependent probability distributions means that the optimal mix 
of spare components at one point in time may not be the optimal mix at another. Thus, 
the approach to take is to compute, for each time period of interest, the marginal in- 
crease in spare parts to achieve a given capability over those already input or determined 
for a previous time. In determining the supply level, the model attempts to provide 
enough spare parts to give the desired confidence at the lowest cost at each point in time 
of interest. Thus, the objective function is to minimize the total cost of spare parts. 

Let: 

S, = the spare parts level for component i 
C, = the unit cost of component i 
a = the desired confidence level 

= the non-mission capable rate not to be exceeded 

P{K^,S) = the probability that the non-mission capable rate is less than K„ 
given a stock level S 

1 = types of repairables required on each aircraft 
Then, the problem to solve is: 



minimize fc,s, 

i=l 



subject to P{K„,S) > a, S/ > S,-o 



where: 

= the input stock level or previous time optimization stock level for component i 
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Assuming complete cannibalization, P{K„,S) equals: 



p(K„,s) = 

/=1 



where: 

P'iQ.K) = Z P.ifO 

K=0 

P,{K) = the probability of exactly K failures of component i 

The necessary’ condition for the performance constraint to be met is to have, 
P‘{QtK„)>a, for each i. Then marginal analysis is used to determine the best mix of 
additional components to achieve the desired goal. This process proceeds by investing 
in one additional component at a time which is selected by finding the component that 
gives the largest increase in the logarithms of the confidence level at the lowest cost. 
That is, we determine: 

A-JnP{K„S)iq 



where: 

MnP{K„,S) = h{P{K„,S‘)l P{K„,S)] 



S' = (S„ S„ .... ,S, .„....) 



The component for which supply is increased one unit is the one whose index solves; 
max A, In P{K„,S) / C,. This process continues until the given confidence level is achieved. 
At this point, the resulting value of S is the efficient solution of the base stockage 
problem. 

(6) Model Evaluation. The Dyna-METRIC model is relatively hard to 
adapt for the Korean Air Force. Especially since this model assumes that each repair- 
able item needing repair immediately receives service. This does not occur in the Korean 
Air Force. Currently, the Korean Air Force applies a batch repair policy for repairables, 
except in case of XORS condition aircraft. This is a serious violation of the assumption 
of Dyna-METRIC. So, the Dyna-METRIC model is difficult to use for the Korean Air 
Force repairable management system without using a continuous repair policy. 
d. Mathematical Models Used In The Navy 

( 1 ) Background. The Uniform Inventory’ Control Program (LTCP) was 
developed in August 1965 to provide a standard system to be used at all U.S Naval 
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Supply Systems Command (XAVSUP) ICPs. The Fleet Material Support OlTice 
(F.MSO) under the direction of XAVSUP is responsible for the system design, ADP 
analysis, programniing and documentation of LTCP. The development of the LTCP 
formulas for inventor}' levels follows the approach used by Hadley and Whitin. [Ref 
13, 19] 

^ 2/ Maihemaiical Assumptions 

1. It is a continuous review system. Wholesale inventory levels requirements and as- 
sets are known by the Inventor}' Control Point (ICP) at all times. [Ref 13, 19; pp. 
162-165, Chapter 3 Appedix A] 

2. It is a steady state environment. The key characteristics of the items managed by 
the ICP are assumed to be constant over the forecast period. These characteristics 
include the forecasted average values and variances of the random variables of the 
rate of customer demand, procurement leadtime, production leadtime, depot repair 
times, depot repair survival rate and the rate of carcass returns. 

3. An order for procurement or repair is placed when the assets reach the reorder level 
or the repair level. 

4. Customer demands and carcass returns do not occur in more than one unit per 
transaction. 

5. The unit procurement cost or repair cost of an item is independent of the magni- 
tude of order quantity or repair quantity. 

6. The cost of a backorder and the time-weighted cost of a backorder can be accu- 
rately quantified. 

7. The reorder level and repair le\'el are always non-negative. 

8. The cost to hold one unit of stock in the inventor}' is proportional to the unit cost 
of the item. 

9. Xo interaction exists among families of items or individual nonfamily items or 
both. Each family or nonfamily item's inventory levels requirements are calculated 
independently of those of families or nonfamily items. 

10. The optimal inventory levels are determined by minimizing an average annual 
variable costs equation composed of order cost plus holding costs plus shortage 
costs. 



UICP Depot Level Repairables (DLRs) Procurement Model. The 
repairables procurement model starts with a total variable cost (TVC) equation which 
is to be minimized. A notable difference from the model for consumables is the inclusion 
of receipt of RFI assets from a repair process in the DLR model. 

(a) Ordering Cost 



Let: 



A = administrative cost to order 
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D = forecasted quarterly recurring demand 
B = forecasted quarterly regenerations 
Q = economic order quantity 



Then, 4(Z) — B)IQ is the expected number of procurements per year. Thus, the annual 
ordering cost is given by 4(Z) — B)AjQ. 



I = percent cost per dollar of inventory held annually 
C = item cost (replacement) 

L = procurement leadtime 
T = repair turnaround time 
R = repair level 

F(x) = the function of probability distribution of leadtime 



The mean demand during a procurement leadtime is given by: 



DT + {D - B) X {L - T) 
Thus, holding cost is given by: 



(b) Holding Cost 



Let: 



demand 




let: 



X = shortage cost per requisition short 
E = military' essentiality weight 
F = quarterly requisition rate 



Then: 



{aE)[A{D - B)IQ]F 
D\ {x - R) F[x) dx 



X>R 
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It equals the cost of the expected number of backorders in a j’ear. The TVC equation 
is svmbolized bv: 



TVC = (D - B)A IQ 

+ IC[_QI2 + R — ((£) — B)xL + BxT} + f (x— /?)F(at) c/x] 
+ (/£)[4(D - B)IQ-] FI Df {x - R)F{x)dx 



Then, by setting dTVC / dQ = 0; 



Q = /[S(Z) - B){A + ;.EF I 



D)l ( 

x>R 



{x - R)F{x) dx I IC] 



Also, by setting dTVC / dR = 0, we find: 




R)F{x) dx = 



QICD 

{QICD + A/.EF{D - 5)} 



Because the expression f ^J^x — R)F{x) dx is the cumulative 
distribution for leadtime demand, this is the quantity defined as RISK. Note that since 
Q and R are related, the reorder quantity Q cannot be solved independently. Thus, 
LTCP approximates Q by using a variation of the economic order quantity formula: 



2 = 




8(Z) - B)A 
IC 



Then, QICD / {QICD + 4?.EF{D — B)} can be computed for RISK determination. 

^4j UICP Depot Level Repairables Repair Model. The repair model also 
starts with a total variable cost equation viewed largely independently of the procure- 
ment problem. Note that its time horizon is the depot level repair turn around time. 
The total variable cost equation for the repair model is the sum of the order cost, holding 
cost and backorder cost. 
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Let: 



Q2 = economic repair quantity 
A2 = repair administrative order cost 
C2 = repair price 
R2 = repair level 

F2{x) = probability distribution of demand during during repair turnaround time 
/-2 = repair shortage cost 

Then: 

TVC = 4 min{D,A)A2 / Q2 

+ 1 C 2 {^ + R2- {D XT) + f {x - R2)F2{x)dx) 

L x> R2 

+ A 2 E {4 I Q2)F I Df (x - R2)F2{x)dx 

x> R2 

By setting 3 TVC / SQ2 = 0: 



Q2 = / 8 min [D,B)A2 + 



,efid\ I 

'^x> R2 



{X - R2)F2{x)dx I IC 2 



Also, by setting 8TVC / = 0 : 

J F 2 {x)dx = Q2IC2D I Q2IC2D + 4 A 2 EFmin{D,B) 

X>0 

Again, Q2 and R2 are related. The UICP model approximates Q2 as follows; 



Q2 = ^/8 min {D,B)A2 / /Cj 

Then, Q2IC2D / Q2IC2D + 4X2EFB can be computed by using the above results. 
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^5; Integrated Repairables Model. As stated earlier, the requirements 
computed by the procurement and repair models are accomplished independently of 
each other. As a result, this leads to a carcass constrained situation. That is, the com- 
puted procurement inventor}' level for an item does not provide sufficient carcasses to 
allow repairs at the computed repair inventor}' level. To solve this problem, Naval 
Supply System Command (NAVSUP) has made some changes. Under the model inte- 
gration, there is only one RISK formula: 

RISK = IC2D I ICjD -b /.FE 

where: 

C3 = (5/DXC2) + (1 - -^){Q 

Also, rather than using a procurement leadtime or a depot level re- 
pair turnaround time, it is uses an average acquisition time as the horizon for computing 
the safety level. The average acquisition time (Lj) is defmed by: 

L 2 = (1 - BID)L + {BID)T 

It is clear that is the weighted average of the procurement leadtime 
and the repair turnaround time because D-B represents the quantity to be procured and 
B represents the quantity from the regeneration. 

f6j Model Evaluation. For the LTCP model, three major problems exist. 
The first problem is that LTCP models is a continuous review inventor}' system. The 
Korean Air Force repairable management system applies a periodic review. The second 
problem is that LTCP assumes that no interaction exists among failures of items. 
Finally, the optimal inventory level for LTCP model is determined by minimizing annual 
variable cost (order costs plus handling costs plus shortage costs). However, the Korean 
Air Force does not have the same objective. 

e. Availability Centered Inventory Model 

^7/ Background. ACI.M was developed by CACI Inc. under the spon- 
sorship of the Ship Support Improvement Project, Naval Sea System Command, 
P.MS-306, and aproved in .March 1981 by the Chief of Naval Operations for use in de- 
termining consumer level stockage quantities for selected equipment [Ref. 20: pp. 1-2]. 
However, this model was initially used as a part of a larger model (Logistics Support 
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Economic Evaluation Model) which provided necessar>^ inputs and enabled comparisons 
with other Na\7 stockage policies. After this, successive refinements for the simplifi- 
cation of the solution procedure and associated computer programs have been made for 
several years. 

This model was originally designed to calculate invent or>' levels for 
all items in the parts breakdown of an equipment and at all stockage facilities in a 
multi-echelon support system. Thus, ACIM is capable of computing levels for all ships, 
intermediate maintenance activities, and depots that use or support the equipment. 
However, this model is mainly used in the provisioning process to compute shipboard 
allowances to achieve specified weapon system readiness levels which have not been 
achieved w'ith the standard protection level models. 

(2j Availablity Measure. ACIM recognizes that the purpose of a supply 
system is to provide sufficient support so that a weapon system is operational when it 
is needed. The terminology used to described this goal is operational availability (A„). 
ACIM defines A^ by the following formula [Ref 21; p. 5j: 

A„ — up time / (up time + down time) 

= MTBF / (MTBF + MTTR + MSRT) 

A^ may also be interpreted as the probability that the equipment is in an operable con- 
dition at a random point in time. Among the three factors of A„, the MTBF and MTTR 
are system parameters outside the control of the supply system and are viewed as con- 
straints. .MSRT is the only term which depends on stockage postures, it is therefore the 
one that the ACLM model focuses on to achieve a given value for A^. 

(3) Mathematical Assumptions 

1. All parts are organized in terms of equipment with a top-dow'n breakdown that can 
be represented as an arborescent network. Any part may be totally consumable, 
totally repairable, or any mix thereof [Ref 21: pp. 10-11] 

2. Stockage and maintenance facilities are organized in a hierarchical structure ac- 
cording to supply/ maintenance flows w'hich can be represented as an arborescent 
network as illustrated below. Each facility has a colocated maintenance and supply 
capability. The facility at the top of the structure is assumed to have an infinite 
supply of all items. 

3. External demands upon supply are stationaiy’ and compound-Poisson distributed. 

4. All stockage locations use a continuous review, (S-1,S) ordering policy. 
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5. Mean Time to Repair (MTTR) is defined to include all equipment downtimes that 
are not supply related. 

6. There is no lateral resupply (transshipment) among bases. 

7. Items repaired at any location are assumed to be returned to collocated stocks for 
issue. 

S. If. for a given facility in the network, the stocks are physically distributed in several 
places, it is assumed that the resupply time for direct customers (next lower ech- 
elon) is independent of the location. 

9. The order policy assumption precludes consideration of economies of scale for re- 
supply. That is, all ordering is on a one for one basis. 

(4) Model formulation. The goal of ACI.M is to maximize the opera- 
tional availability {A^) of a weapon system subject to a given inventory budget. With the 
definition given above of A, and with the view that MSRT is the onlv term which is af- 
fected by the stockage decision, the developers of ACI.M argue that the allocation which 
maximizes A^ is equivalent to the allocation which minimizes MSRT. ACIM therefore 
actually attempts to minimize .MSRT subject to given constraints. 

Let i be an arbitrarx’ item in equipment e (which may be e itselQ- Let 
u = 0 represent an arbitraiy facility in the support system and u= 1,2,3,...., represent fa- 
cilities at the next lower level, i.e. those facilities that submit items for repair directly to 
or obtain resupply from facility o. Then, the objective can be explicitly stated as follows: 

Find values for which minimize for all user locations u: 



subject to 

k,v 

where: 

C* = unit cost of item k 
B = given budget for spares procurement 

D,^ = expected delay per demand upon inventory’ for for item i at location u 
5*, = the level of inventory for item k at location v 

D,^ is one of the components of which is defined as the mean time to return a failed 
unit of item i at location u to a serviceable condition. is given by: 

Mju = Z),„ + Tju 
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where: 

= expected delay per demand upon inventory for for item i at location u 
= mean time to repair item i at user location u 

Z),„ in turn is given by: 



A„ - (I / (■'' - AJ P (X ■. 

where: 

S,„ = stock level of item i at location u 

= expected number of demands upon inventory for item i at location u 
P{X : /,„rj = probability of X units of stock reduction for item i at location u 
r,„ = mean resupply time for item i at location u 

T;„ is given by: 

T^iu = i^u + Uu) + ( 1 ~ l'iu)i^iu + ^iu) 



where: 

= probability that a demand for item i upon inventory at location u 
L,„ = average resupply leadtime assuming stock is is available 
Ll = additional resupply leadtime due to expected shortage at the resupply source 
= average shop repair cycle assuming availability of spares for items within i 
= additional shop repair cycle due to expected shortage of spares items within i 

Since ACIM assumes that MTBF and MTTR are independent of the 
stockage policy, minimizing is assumed to maximize (operational availability of 
a equipment at location u). 

(5J Solution Technique. Assume that initial values for S,„ are given for 
all items and locations. These may all be zero or some minimum value given by policy 
or current assets. Then compute the MSRTs for all items. 

ACIM calculates what the new MSRT would be if one additional 
unit of stock were placed against that particular item. Subtracting the new MSRT from 
the old MSRT and then dividing by the unit cost of the item provides the model with a 
value which is multiplied by the item's demand to determine a selection-rank. After the 
selection number is calculated for each item, one unit of stock is added to the candidate 
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with the highest selection rank number. The operational availability (AJ is calculated 
and if the target A^ has not been achieved, the model continues the same procedure. 

Model Evaluation. As with the LTCP model, three major violations 
of the basic assumptions of the ACl.VI model exist. The first is that ACIM assumes that 
you have a continuous review system. The second is that ACIM assumes there is no 
lateral resupply. Finally, ACIM assumes that the stocks are physically distributed at 
several locations. Obviously, the Korean Air Force supply system does not require se- 
veral physical locations, even though the system objective is to maximize operational 
availability. 

D. SUMMARY 

Throughout Section B and C of this chapter, we have reviewed the various inven- 
tor>’ models and theories for repairable items which are currently in use in the inilitary 
services, or in the process of being evaluated. All of the assumptions and the math- 
ematical formulations of the models were stated and discussed. 
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IV. MODEL PROPOSAL FOR THE KOREAN AIR FORCE 



A. INTRODUCTION 

In Chapter II, we summarized the current Korean Air Force repairables manage- 
ment system and repair process. The general structure of repairables management in the 
Korean Air Force is similar to the structure assumed by several of the models described 
in Chapter III. 

In this Chapter, we describe a proposed inventory control process for the Korean 
Air Force. This process is based partly on queueing models, and partly the Wilson- 
Harris economic order quantity model. 

B. THEORETICAL BASE FOR THE MODEL PROPOSAL 

1. Notation 

c = number of identical repair channels 

/X, = average repair rate per repair channel 

a, = average failure rate of an individual unit of a repairable item 

iV, = random variable describing the steady 
state number of XRFI units of item i 

P„ = steady state probability that there are n NRFl units 

2. The M/M/l/K/K Queueing System 

This model is a limited source queueing model in which there are only K cus- 
tomers. It is variously called the machine repair model, the machine interference model, 
or the cyclic queueing model. It is one of the most useful of all queueing theory models. 
One way to view this model is shown in Figure 6. [Ref 22: pp. 186-190] 

The population of potential customers for this queueing system consists of K 
identical devices, each of which has an operating time of O time units between break- 
downs, O having an exponential distribution w'ith average value 1 / a. The one 
repairman repairs the machines at an exponential rate with an average repair time of 
1 / pL time units. The operating machines are outside the queueing system (outlined by 
the dashed lines) and enter the system only when they break down. The queueing system 
always reaches a steady state because there can be no more than K customers in the 
system (one machine being repaired and K-1 waiting for repairs). 
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Figure 6. Queueing System nitli One Keuainnan (M/M/J/K/IQ 



NMicn n oftlic macliincs arc down (not operating), tlicn K-n of tlicin arc oper- 
ating and the time until tlie next niacliine breaks down is the minimum of K-n identical 
exponential distributions with parameter (K-n) a. 

3. The M/M/c/K/K/ Queueing System 

This queueing system is siinilar to the machine repair model considered in the 
section above except that we have e rather than one repairman, where c < K, tis shown 
in Figure 7. [Ref. 22: pp. I90-192J 

C. DEVELOPING A MODEL FOR THE KOREAN AIR FORCE 
1. Application Conditions 

1. Supporting I) items, i.c.. we want D items to be available at all times. In this ease. 
D represents the number of aircraft which the Korean Air Force wants to have 
available at all times. 

2. P(Succcssful Repair) is close to 1. When P(succcssrul Repair) = 1, we have a pure 
machine repair queueing system. When P(Sucecs.srul Repair) = 0. we have a eon- 
sumablc item whose management is best accomplished using standard methodol- 
ogies from inventory theory. The sample ilata collected from the Korean /\ir Force 
for this study indicated that P(Succcs.sful Rc]xrir) = 0.9S or higher for all rcjxiirable 
items examined. For this reason, the proposed model is based partly on the m.i- 
chinc repair qucueiiig model. The data base used for the study was obtained from 
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Figure 7. Queueing System wilii c Repnirman (M/iM/c/K/K) 



the KAFLC during the 19S4-19SS. Table 3 shows tlie mean probability of sue- 
cesslul repair based on the sample dataA 



Table 3. P(SUCCESSFUL REPAIR) COMPUTED FROM SAMPLE DATA 



Year 

A'C 


S4 


85 


86 


87 


88 


Average 


F-16 


NA 


NA 


1 .000 


I.Onn 


1.000 


1.000 


F-5 


U.9S8 


0.980 


0.987 


0.986 


0.984 




F-4 


0.986 


0.987 


0.986 


0.976 


0.987 


0.984 



3. Conditions close to the M/M/l. K/K queueing system hold. Because, there are ;i 
number of repairmen in repair depot, it might be reasonable at First to assume tluit 
the M/M/c;K/K mudcl applies. However, most repairs require spceiali/cd skills 
held by few people for a given repairable part. Further, the rciiair depot is stalled 
and scheduled so that there is little overall excess capacity. As a result, there is 
little change in the repair process output rate when the second, third, fourth, etc. 
unit backs up in the repair queue. 

4. Steady state has been achieved. 



6 The sample data list appears in Appendix A and B 
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2. Model Formulation 

This proposed model consists of 3 steps. The first step computes a minimum 
population size A' for the iih type of repairable item. Population refers in this case to 
the installed material plus the stored RFI material plus the NRFI material awaiting re- 
pair. In other words, all material, both good and broken. 

Once the minimum population size K, is chosen, an attrition buy safety level 
S, is chosen. This safety level is used to insure that the actual population size usually 
stays above A' during the period when the system is awaiting the arrival of an attrition 
buy. 

The last step is to compute the size of an attrition buy itself. This Q', is com- 
puted using the Wilson-Harris economic order quantity model. 

In step 1, described on the next page, the M .M 1,K,K queueing system will 
be used to obtain probability e.xpressions of the form P(.V, > A' — Z),). These expressions 
are obtained as follows: 



K, 

P{\\>K--D^= ^ P{h\ = n) 

n- - D, +1 



where: 

.V, = Discrete random variable representing the number 
of units of item i which are broken 

A' = Decision variable representing planned minimum 

number of units of the item to be owned, including those 
installed in the repairable equipment 

D, = Number of units of item i which the Korean Air 
Force wants to have available at all times 

The right hand side of the last expression is equal to: 

K,-D, 

1 - ^ = «) 

«=o 

using the M'MT,K/K queue, Allen [Ref 22: p. 192] gives the following formula for 



P(.V, = /r) 



Af ! a, _ 
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where: 



= lE 




k=0 

In step 2, described in the following section, the expression P{Xj> S,) is also 
used. It is obtained using the Poisson distribution: 



A' = discrete random variable representing attrition 

demand for item i during a procurement leadtime 

D^i = forecasted attrition demand for item i 

Li = forecasted procurement leadtime for item i 

We will want S, such that P(A’, > 5,) is less than a risk value, b 2 „ which is pre-specified 
by the user, where P{Xi > S,), is given by: 



P(A; = A-) = 




where: 



oo 






k=S,+\ 




oo 



k= Si +1 




5 , 



a. Step I 

There are several different objective functions which might be used to help 
choose the size of the population fot the ith item. The following four subsections de- 
scribe four different math programs for choosing Aj. 

(]) Alternative I 



Choose A' by the following: 
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min Ki 

subject to P(.V;- > (A"; — Dj)) < by, 



where: 

D, = Number of units of item i which the Korean Air 
Force wants to have available at all times 

by = Parameter chosen by the decision maker 

This first alternative formulation contains no explicit budget constraint. However, 
minimizing K, for all items is equivalent to minimizing the cost of that portion of the 
population which is owned to meet the availability goals. Unfortunately, this formu- 
lation doesn't identify solutions which exceed the available budget. The next math 
programming formulation accomplishes this. 

^2) Alternative II 



m 

min 

i=l 

subject to Ki — Di)<by, /=1, ,ni 

m m 

and ^^hjKi -f ^<3,[max(0, A; — Ao,-)] < B 

i=\ i=l 

where: 

h, = annual holding cost per unit held 

B = total budget for one year 

m = number of repairable items being considered 

A'„ = previous planned minimum number of units owned 

a, = unit cost of item i 

Alternative III 

Maximize customer service subject to investment constraint: 
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m 



min Z = ^P(A>(A^ - A)) 
/=1 



m m 

subject to '^^hiKi + ^fl,[max(0, Ki — < B 

i =\ /=1 

(4) Alternative IV 

Maximize customer service subject to investment constraint (another form): 



min Z = max {P(A)>(A,- - Z);))} 

m m 

subject to y^hjKi + ^a,[max(0, A) — ^ 5 

(=1 j=i 

The choice between these alternate formulations depends on the way the Korean Air 
Force feels most comfortable with resolving the conflict between budget limitations and 
customer service. For the numerical computation of this study, we chose alternative 11. 
This choice balanced the need for realism against the time constraints imposed on this 
research. 

b. Step 2 

Attrition buy reorder point, A, = A', + S„ where S, is obtained from the 

following; 



min Si 

subject to P{Xi > Si) < b 2 i 



where: 

A] = discrete random variable representing attrition 
demand for item i during a procurement leadtime 

^ 2 , = parameter chosen by the decision maker 
for item i 

Using the reorder point A„ place an attrition buy when the total number of unattrited 
units (working and failed, installed and uninstalled) reaches or drops below A,. 
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c. Step 3 

Use the Wilson-Harris EOQ model to choose procurement quantity: 



At = administration and ordering cost for item i 
D., = attrition demand rate for item i 
It = annual holding cost rate for item i 
C, = unit cost for item i 

So the initial buy quantity formula is + Q', = K, + S, + Q‘. 



This proposed repairable item inventory control process works as follows. 



To start the system, purchase A' + S, + O', units. As installed items fails, repair them 
as soon as possible (i.e. do not use batch repair). When attritions in the repair process 
reduce the actual population size to K, + S„ place an attrition buy. 

D. SUMMARY 

In this chapter, we developed a new inventory model based on the M .MT/K/K 
queueing theory and the current repairables management system of the Korean Air 
Force. In the next chapter, we will compare the current inventory model with the pro- 
posed model, in terms of the numerical results of each, using computer computations. 




where: 
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V. COMPARING THE COMPUTATION RESULTS OF THE MODEL 



A. INTRODUCTION 

In this chapter, we will discuss a specific result for the two models, one is a existing 
model of the Korean Air Force, the other is a proposed model. The results show only 
summary of the computer computation and a comparison of the results. The whole 
computation results are shown in Appendix C. 

The purpose of this chapter is to provide a numerical analysis between the proposed 
model and the current Korean Air Force model. The data used in the computer com- 
putation consists of a selection of 8 line items form the data on 50 line items which was 
obtained from the Korean Air Force. Due to the limited time available for this research, 
only 8 items were selected for use in the comparison of the models. A selection of high 
demand and low demand items were chosen. The complete data set is described in Sec- 
tion C of this chapter. 

B. MEASURE OF EFFECTIVENESS (MOE) 

The following .MOEs are all reasonable methodologies for evaluating repairable item 
inventory models for the Korean Air Force. How'ever, we will concentrate on one MOE 
for model comparison such as operational availability. 

1. Fill Rate 

It is computed by taking the total number of units demanded at a base over a 
fixed period time and dividing that number into the total number of units issued at the 
time they were are demanded. Thus, fill rate is the percentage of demands that are filled 
immediately. 

2. Operational Rate 

Operational rate is the probability that, at any given point in time, there will be 
no stockouts from base supply (backorders). Operational rate is computed by counting 
up the length of time (in days) during a year that no backorders existed and dividing this 
number by 365. This gives us the percentage of time during the year that no backorders 
were in existence. 

Operational rate has an advantage over both the fill rate and mean number of 
backorders in that it may be directly related to the supply system's effect on operations. 
However, it has a disadvantage over both the fill rate and mean backorders in that it has 
a rather bothersome all-or-nothing character. 
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3. Mean Supply Response Time (MSRT) 

MSRT is the mean time it takes for the supply system to respond to the demand 
for a replacement part or component. It is obtained by taking each stockout that is es- 
tablished during a period of time, observing how many days it takes to satisfy the 
backorders, adding up all these numbers and dividing the sum by the total number of 
demands during that period. It is the average time weighted unit short (TWUS). 

.MSRT is important by itself as an indicator of success of the supply system in 
meeting response time goals. As a measure, it considers both the likelihood of satisfying 
demands from stock on hand and the length of the delay in satisfying demands when the 
system runs out of stock. 

-1. The Average Number of NORS Aircraft 

NORS stands for "not operationally ready because of supply". The average 
number of NORS aircraft is computed by counting for each aircraft the number of 
NORS days during the course of a year, adding these numbers up for each aircraft, and 
finally dividing by 365. Considering that the purpose of a spare parts supply system for 
the Korean Air Force is to maintain the operational readiness of aircraft, the average 
number of NORS aircraft would certainly seem to be a reasonable measure of efiective- 
ness. A stockage model attempting to optimize with respect to a NORS measure would 
require more restrictive assumptions than those optimizing with respect to fill rate, the 
average backorders, MSRT. or the operational rate. Above all, its lack of "separability" 
is a cause of major mathematical problems. 

Something to note is that the operational sector of the KAF has consistently 
suggested that the Korean Air Force supply system should be managed in operational 
terms, not by supply terms. That is, their imminent question has been "How many 
spares are needed to provide an operational readiness of X percent of the aircraft?". 

5. Constraints 

The budget available for the investment in spares is either a constraint in all of 
the models or the objective is to achieve specified performance at the minimum budget. 
The essence of any inventor}' control problem is the trade-off between cost and system 
performance. As mentioned in the beginning of Chapter II, the budget constraint is the 
major resource constraint for the Korean Air Force. 
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C. COMPARING THE MODEL RESULTS 



1. Data Characteristic 

The data used for the model comparison are 8 line items of the 50 line items in 
the sample data. The whole sample data is shown in Appendices A and B. This data 
consists of; a) quarterly demand data for 19 quarters; and b) quarterly attrition data for 
19 quarters for each item. The demand data was analyzed using Minitab. The mean, 
standard deviation, variance and squared coefficient of variation were calculated for each 
item. Since data concerning the failure process for each installed item wasn't available, 
we estimated the mean item failure rate, a, by computing mean demand over 5 years, 
divided by number of aircraft. We estimated the annual item repair rate, /r„ by dividing 
actual working hours per year by the average repair time in hours for the item. The 
current average repair time for each of the 8 items in the sample was obtained from the 
Korean Air Force [Ref. 23]. These values are shown in the table below. For the pur- 
poses of this comparison we assumed that the random variables associated with the 
failure process and the repair process had exponential distributions. 



Table 4. MEAN REPAIR TIME: 



This table shows the mean 
active repair times for the 8 
sample items 



Item number 


Mean Repair 
Time (hours) 


1 


2.9 


2 


5.3 




3.9 


4 


2.1 


5 


7.5 


6 


6.1 


7 


1.5 


8 


1.2 



The actual data for an item consists of 19 quarterly demand observations during 
1984-1988. Items were selected from three different kinds of aircraft, the F-4D/E, 
F-5A/B and F-16A/B. Appendix B shows the Minitab results for all 50 items in the 
sample data set. 
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The data used in this research also included quarterly attrition quantities for 19 quarters. 
However, for the model comparison, we selected 8 items." The items which were chosen 
had a variety of annual demands. 

2. Comparing the Results 

The complete computation results for both models are shown in Appendix C. 
A summary of these results is shown on the following pages. For the Korean Air Force 
model, the target operational availability is 90 percent. In the Korean Air Force, this 
is a fixed level of availability. 

During the computational steps, the actual operational availability allows an 
increase from 90 percent to 93 percent, however the proposed model didn't meet the 
budget constraint for a 94 percent operational availability. Figure 8 shows the difference 
of K values; a comparison of the current K values of the Korean Air Force with min K 
from the proposed model. Figure 9 shows the difference in safety levels for the proposed 
model and the Korean Air Force model. In the comparison, the attrition protection 
levels of 0.98 is fixed. Appendix C shows the results based on several attrition protection 
levels between 0.90 and 0.98. SLQ is a safety level for the Korean Air Force model, min 
S is the attrition safety level for the proposed model. Figure 10 shows the difference in 
procurement quantities. For the Korean Air Force model, RO is a procurement level, 
optimal Q is for the proposed model. In the actual difference for the two models, the 
proposed model has slightly lower procurement quantities than the Korean Air Force 
model. 

3. Description of the Computer Program 

For the Korean Air Force repairable item inventor}' model, we first computed 
a demand forecast using a trend value to guide our choice of alpha value to be used in 
an exponential smoothing forecast for the next quarter's demand. When the trend value 
fell between 0.9 and 1.1, we chose an alpha value of 0.2, otherwise we used an alpha of 
0.4. Then the demand forecast was used to calculate the DDR, DRP, RCQ, OSTQ, SLQ 
and RO. The computational procedure for this current Korean Air Force model is 
briefly described in the following paragraph. The complete computational results for the 
Korean Air Force model are shown in Figures 20 to 27, and the whole mathematical 

~ hem numbers for the 8 items chosen from the 50 items are 10, 16, 20, 27, 30, 33, 37 and 
46. 

8 Attrition protection level is the probability that attrition will not reduce the population of 
item i below K, when an attrition reorder is placed when the population drops to the A', + 5, level. 
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Figure 8. Difference of Total Stock Level 
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Figure 9. Difference in Safely Levels Det\\eeii the Two Models 
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Figure 10. Dirfcrencc of Optimal Procurement Quantity 
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process for the Korean Air Force repairable item inventor}' model is described in Chap- 
ter 11. 

For the DDR, we computed the total demand for each quarter, divided by 365 
days. The DRP equals the RTS divided by the sum of RTS, NRTS and COND. The 
RCQ was computed by multiplying DDR times DRP times RCT. This RCT refers to 
the repair cycle time allowed for depot level maintenance. For these computations, we 
used 120 days for RCT because this is the ma.ximum allowed by the Korean Air Force. 
The OSTQ is computed by multiplying DDR times NDRP times OST. For OST, we 
used 365 days which is the maximum allowed by the Korean Air Force. The SLQ is 
computed by finding the square root of 3 times RCQ plus OSTQ. The RO is computed 
by adding SLQ, RCQ and OSTQ. 

For the model comparison, we used two system parameters (SLQ and RO) and 
the item's current K value. Using the historical demand data obtain from the Korean 
,Air Force, and their current repairable item inventor}' as described above, the requisition 
objective (RO) and safety level quantity (SLQ) were computed for each calendar quarter 
in 19SS (note that the equivalent historical information was not available from the 
Korean Air Force). The current K value in Figure S was obtained by taking the sum 
of RFl, NRFl and the currently installed population, as of November 19SS (previous 
historical data was unavailable). The SLQ for model comparison in Figure 9 is the av- 
erage safety level in 19SS. Also, the RO in Figure 10 is the average RO level in 19SS. 

For the proposed model, we first need to choose for the first constraint set 
for each item. For by, one minus the current target operational availability in the 
Korean Air Force was used. To obtain the constraining total budget value B, we com- 
puted the sum of the current K,, value for each item times its unit cost plus the holding 
cost for each item, over all 8 items. The rho value ( -rp ) for an item was derived from 
the average failure rate (x) of that repairable item, divided by its average repair rate (^). 

For the model comparison, a value for the budget constraint had to be obtained. 
We estimated a value for this variable by summing the following products over all items: 
K,! times unit cost plus annual holding cost for item i (reminder: K^, is the number of 
units of item i currently owned by the Korean Air Force). 

In step 1, the min K, values are computed based on two constraints. First, we 
calculate the probability that .V, is greater than K, — D, for an item. The starting value 
used for A', is D, (the size of the currently installed population). For the min A, in the 
computer program, we first compute of the M/M/l/K/K queueing system. And then, 
we compute P{S\ = «) using P^. In the P{\\ = n) computation, n equals K, — D,. If the 
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constraint on P{i\\ > K, — D,} isn't satisfied by choosing K, = Z), then vve increased K, by 
1. These steps are continued, until /*(A’ > A' — D,} is less than 0.10. At this time, we 
took the current n value, and then added to get the min A', value. After obtaining the 
min A' values from the first constraint set for each of the 8 items, the budget constraint 
is checked. If these A' values don't satisfy the budget constraint, we go back to the first 
constraint set and increase the b^| values. In this case, a method for choosing an in- 
creased set of by 's must be chosen. Two ways exist to do this. The first method is to 
increase the by values simultaneously for all values. The other way of solving this 
problem is to apply different operational availabilities for each item. In other words, the 
operational availability chosen would depend on the item characteristics. The first 
method described above (simultaneous increase of ever\' by) was used in this work. If 
the min satisfies the budget constraint, this is min A’] used the model comparison, in 
Figure 8. As it turns out, the first set of min A, values met the budget constraint. So, 
the min K, values from the first constraint set gave a feasible solution in the budget 
constraint. Actually, the budget constraint allowed up to 93 percent operational avail- 
ability for all 8 items. 

For step 2. we use the Poisson probability distribution to get minimum attrition 
safety levels. Suppose that risk ^ 2 > = 0.1. This would mean that the attrition protection 
level is 0.90. Therefore, 1 — ^ 2 , equals the protection level. For the min S„ we first cal- 
culate probability of A”, = 0, and next calculate probability X, = 0 and so on. We sum 
these individual probabilities up until the cumulative value meets or exceeds the pro- 
tection level (for the model comparison, we use 0.98). This is the min value for the 
model comparison in Figure 9. 

For step 3, the Wilson-Harris EOQ formula is used. For this formula, a per 
order ordering cost of 10 percent of the unit cost was used. This is the value which the 
Korean Air Force currently believes to be in effect [Ref 1: p. 65). Similarly, a holding 
cost rate of 10 percent per year of unit cost was used. The resulting optimal procure- 
ment quantities for the model comparison are shown Figure 10. The complete compu- 
tational results for the proposed model are shown in Figures 28 to 35. 

4. Suniinar}' of the Results 

The proposed model allows for a 93 percent operational availability within the 
budget constraint. This is an improvement of 3 percent in operational availability with 
the same budget that was used with the current Korean Air Force model to obtain 90 
percent operational availability. 
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D. SUMMARY 

In this chapter, we compared two models based on computer computations using a 
sample data set. The proposed model has shown an improvement over the current 
model. However, it doesn't mean the proposed model is better for the Korean Air Force. 
A tremendous efibrt would be required to actually implement this model in the Korean 
Air Force. Furthermore, additional research should be done to explore the performance 
of the proposed model. 
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VI. SUMMARY AND CONCLUSION 



For this research, there were three major objectives. The first was to review the lit- 
erature to understand the current Korean Air Force model and other models, such as 
METRIC, Mod-METRIC, for controlling the stockage decisions in multi-item, multi- 
echelon inventory systems involving repairable items. The second objective was to de- 
velop a new model for the Korean Air Force inventory management system. The last 
objective was to evaluate this new model by comparing the current Korean Air Force 
model with this proposed model. 

For the first objective, we studied the METRIC family model, ACIM model and 
U.S Na\w model. We then reviewed their basic assumptions and mathematical ap- 
proach. For the second objective, we reviewed the Korean Air Force inventory man- 
agement system concept and conditions, and developed a proposed model. In order to 
accomplish the last objective, sample data from the KAF was used in comparing the 
Korean Air Force model and proposed model. 

Analysis of the sample data revealed that the proposed model yielded some im- 
provement in operational availability and safety level with the same budget. As men- 
tioned in Chapter V, the proposed model requires further research. Currently the 
research is hard to do because of time constraints and data available from the Korean 
.■Mr Force. So, additional research should be done to eliminate these shortcomings in the 
proposed model. Also, for the implementation, we need more data points to obtain 
better distributions for the demand process and repair time process. 
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APPENDIX A. SAMPLE DATA LIST 1 



A. SAMPLE DATA DEMAND DISTRIBUTION LIST 



Table 5. DEMAND DISTRIBUTION LIST 



QTR 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


84-1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


84-2 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


84-3 


0 


0 


0 


0 


0 


0 


0 


0 


0 


2 


84-4 


0 


0 


0 


0 


0 


0 


0 


0 


0 


2 


1984 

Total 


0 


0 


0 


0 


0 


0 


0 


0 


0 


4 


85-1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


1 


85-2 


0 


0 


0 


0 


0 


0 


0 


0 


0 


2 


85-3 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


85-4 


0 


0 


0 


0 


0 


0 


0 


0 


0 


2 


1985 

Total 


0 


0 


0 


0 


0 


0 


0 


0 


0 


5 


86-1 


0 


0 


0 


2 


0 


0 


0 


0 


0 


0 


86-2 


0 


0 


0 


0 


0 


0 


0 


0 


0 


2 


86-3 


0 


0 


0 


1 


0 


1 


0 


0 


1 


0 


86-4 


0 


0 


0 


0 


0 


1 


0 


0 


2 


3 


1986 

Total 


0 


0 


0 


3 


0 


2 


0 


0 


3 


5 


87-1 


0 


2 


0 


2 


1 


0 


0 


0 


5 


2 


87-2 


2 


1 


0 


4 


0 


1 


1 


1 


7 


0 


87-3 


0 


0 


0 


0 


0 


0 


0 


1 


2 


0 


87-4 


1 


0 


1 




0 


2 


0 


1 


2 


4 


1987 

Total 


3 


3 


1 


9 


1 


3 


1 


3 


16 


6 


88-1 


0 


2 


1 


1 


0 


0 


2 


0 


4 


3 


88-2 


0 


1 


0 


3 


0 


3 


0 


1 


6 


0 


88-3 


1 


0 


0 


5 


1 


0 


0 


1 


4 


2 


1988 

Total 


1 


3 


1 


9 


1 


3 


2 


2 


14 


5 
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Table 6. DEMAND DISTRIBUTION LIST 



QTR 


11 


12 


13 


14 


15 


16 


17 


18 


19 


20 


84-1 


0 


0 


0 


1 


0 


6 


1 


2 


1 


0 


84-2 


0 


0 


0 


2 


1 


1 


0 


4 


7 


3 


84-3 


0 


0 


0 


2 


0 


0 


0 


0 


2 


0 


84-4 


0 


0 


0 


2 


0 


2 


0 


2 


5 


2 


1984 

Total 


0 


0 


0 


7 


1 


9 


1 


8 


15 


5 


85-1 


1 


0 


4 


0 


1 


1 


0 


8 


4 


2 


85-2 


0 


0 


0 


5 


0 


1 


0 


0 


3 


0 


85-3 


0 


0 


0 


2 


0 


3 


0 


0 


1 


2 


85-4 


0 


1 


0 


3 


1 


0 


0 


0 


1 


3 


1985 

Total 


1 


1 


4 


10 


2 


5 


0 


8 


9 


7 


S6-1 


0 


1 


0 


4 


0 


0 


2 


2 


3 


4 


86-2 


0 


0 


0 


1 


1 


2 


0 


5 


3 


0 


86-3 


0 


0 


0 


3 


3 


0 


0 


1 


5 


3 


86-4 


0 


0 


0 


2 


1 


1 


0 


1 


1 


2 


1986 

Total 


0 


1 


0 


10 


5 


3 


2 


9 


12 


9 


87-1 


4 


0 


0 


5 


2 


1 


0 


1 


0 


0 


87-2 


1 


1 


0 


4 


5 


0 


0 


2 


4 


0 


87-3 


5 


0 


1 


1 


3 


1 


0 


0 


1 


1 


87-4 


2 


0 


0 


'y 


2 


1 


0 


0 


2 


2 


1987 

Total 


12 


1 


1 


13 


12 


3 


0 


3 


7 


3 


88-1 


1 


0 


1 


4 


3 


0 


1 


0 


2 


0 


88-2 


2 


1 


0 


2 


2 


1 


0 


1 


3 


0 


88-3 


1 


0 


0 


4 


0 


0 


0 


0 


1 


2 


1988 

Total 


4 


1 


1 


10 


5 


1 


1 


1 


6 


2 
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Table 7. DEMAND DISTRIBUTION LIST 



QTR 


21 


22 


23 


24 


25 


26 


27 


28 


29 


30 


S4-1 


0 


2 


0 


1 


9 


4 


36 


84 


4 


4 


84-2 


1 


0 


0 


1 


5 


5 


44 


66 


7 


0 


84-3 


1 


0 


0 


0 


11 


4 


52 


43 


1 


0 


84-4 




2 


4 


3 


1 


0 


55 


32 


7 


3 


1984 

Total 


5 


4 


4 


5 


26 


13 


187 


225 


19 


7 


85-1 


1 


4 


1 


10 


0 


2 


101 


71 


2 


3 


85-2 


1 


1 


0 


0 


10 


2 


67 


59 


0 


0 


85-3 


1 


1 


0 


5 


5 


5 


28 


50 


5 


0 


85-4 


0 


1 


1 


4 


4 


1 


42 


25 


2 


2 


1985 

Total 


3 


7 


2 


19 


19 


10 


238 


205 


9 


5 


86-1 


2 


5 


4 


5 


7 


6 


61 


40 


0 


3 


86-2 


4 


3 




1 


8 


1 


17 


39 


4 


1 


86-3 


0 


2 


1 


1 


8 


1 


95 


12 


5 


0 


86-4 


1 


4 


0 


4 


4 


3 


47 


73 


0 


5 


1986 

Total 


7 


14 


8 


11 


27 


11 


220 


164 


9 


9 


87-1 


5 


6 


2 


2 


10 


2 


19 


50 


11 


' 3 


87-2 


4 


1 


0 


2 


7 


7 


24 


37 


7 


0 


87-3 


1 


3 


0 


0 


3 


0 


25 


21 


4 


0 


87-4 


0 


1 


7 


0 


0 


4 


70 


24 


13 


4 


1987 

Total 


10 


11 


9 


4 


20 


13 


138 


132 


35 


7 


88-1 


0 


5 


1 


5 


9 


10 


42 


31 


5 


2 


88-2 


6 


2 


0 


6 


9 


1 


21 


18 


9 


0 


88-3 


0 


3 


2 


2 


5 


3 


18 


23 


3 


4 


1988 

Total 


6 


10 


3 


13 


23 


14 


81 


72 


17 


6 
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Table 8. DEMAND DISTRIBUTION LIST 



QTR 


31 


32 


33 


34 


35 


36 


37 


38 


39 


40 


84-1 


0 


0 


0 


9 


1 


2 


2 


0 


12 


1 


84-2 


5 


0 


1 


7 


7 


1 


3 


0 


21 


7 


84-3 


4 


0 


2 


9 


0 


1 


0 


0 


6 


3 


84-4 


0 


0 


1 


18 


2 


9 


3 


5 


2 


1 


1984 

Total 


9 


0 


4 


43 


10 


13 


8 


5 


41 


12 


85-1 


12 


0 


1 


17 


10 


4 


2 


1 


17 


2 


85-2 


20 


3 


1 


2 


26 


3 


0 


1 


9 


0 


85-3 


1 


0 


3 


2 


32 


1 


1 


0 


13 


5 


85-4 


14 


0 


0 


7 


25 


3 


2 


2 


5 


1 


1985 

Total 


47 




5 


28 


93 


11 


5 


4 


44 


8 


86-1 


17 


1 


0 


4 


30 


5 


2 


4 


10 


0 


86-2 


25 


0 


1 


2 


11 


0 


5 


0 


0 


13 


86-3 


0 


1 


0 


9 


17 


0 


5 


0 


5 


7 


86-4 


21 


2 


1 


0 


4 


6 


0 


3 


6 


6 


1986 

Total 


63 


4 


2 


15 


62 


11 


12 


7 


21 


26 


87-1 


5 


0 


0 


4 


9 


3 


1 


1 


2 


5 


87-2 


5 


1 


2 


5 


5 


3 


2 


1 


2 


5 


87-3 


0 


0 


2 


5 


2 


1 


0 


0 


7 


5 


87-4 


4 


0 


0 


0 


17 


7 


3 


0 


4 


3 


1987 

Total 


14 


1 


4 


14 


33 


14 


6 


2 


14 


18 


88-1 


10 


0 


0 


0 


0 


1 


0 


0 


0 


0 


88-2 


3 


0 


1 


0 


10 


1 


0 


0 


4 


3 


88-3 


7 


0 


1 


5 


4 


2 


4 


0 


5 


0 


1988 

Total 


20 


0 


2 


5 


14 


4 


4 


0 


9 


3 
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Table 9. DEMAND DISTRIBUTION LIST 



QTR 


41 


42 


43 


44 


45 


46 


47 


48 


49 


50 


S4-I 


0 


0 


1 


4 


7 


14 


28 


5 


3 


13 


S4-2 


0 


0 


0 


9 


7 


0 


30 


6 


3 


9 


S4-3 


0 


0 


0 


9 


9 


0 


17 


11 


9 


10 


S4-4 


0 


0 


0 


10 


2 


11 


29 


12 




17 


19S4 

Total 


0 


0 


1 


32 


25 


25 


104 


34 


18 


49 


85-1 


8 


1 


2 


5 


1 


5 


23 


8 


5 


15 


85-2 


8 


1 


2 


5 


1 


10 


23 


8 


5 


15 


85-3 


1 


7 


0 


8 


4 


0 


13 


13 


0 


8 


85-4 


3 


3 


4 


5 


2 


6 


11 


7 


4 


8 


1985 

Total 


20 


12 


8 


23 


8 


21 


70 


36 


14 


46 


86-1 


4 


0 


0 


10 


0 


9 


7 


12 


2 


20 


86-2 


7 


4 


1 


13 


11 


0 


19 


10 


3 


0 


86-3 


1 


1 


1 


7 


9 


0 


18 


5 


7 


6 


86-4 


4 


2 




5 


5 


11 


6 


16 


0 


16 


1986 

Total 


16 


7 


5 


35 


25 


20 


50 


43 


12 


42 


87-1 


0 


5 


0 


1 


6 


9 


5 


9 


0 


8 


87-2 


0 


1 


1 


1 


6 


7 


9 


10 


6 


1 


8"-3 


0 


1 


0 


9 


2 


0 


8 


3 




9 


87-4 


2 


0 


0 


6 


4 


13 


8 


11 


0 


15 


1987 

Total 


2 


7 


1 


17 


IS 


29 


30 


33 


9 


33 


88-1 


0 


0 


0 


0 


5 


0 


7 


11 


0 


4 


88-2 


0 


4 


1 


3 


2 


0 


5 


4 


5 


7 


88-3 


1 


0 


0 


3 


1 


12 


14 


7 


1 


6 


1986 

Total 


1 


4 


1 


6 


8 


12 


26 


22 


6 


17 
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APPENDIX B. SAMPLE DATA LIST 2 



A. SAMPLE DATA DISPOSAL RATE DISTRIBUTION LIST 



Table 10. DISPOSAL RATE DISTRIBUTION LIST 



QTR 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


84-1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


84-2 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


84-3 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


84-4 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


85-1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


85-2 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


85-3 


0 


0 


0 


0 


0 


0 


0 


0 


0 


1 


85-4 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


86-1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


86-2 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


86-3 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


86-4 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


87-1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


87-2 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


87-3 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


87-4 


0 


0 


0 


0 


0 


0 


0 


0 


0 


1 


88-1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


88-2 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


88-3 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


P(successfu 

repair) 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


0.92 
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Table 11. DISPOSAL RATE DISTRIBUTION LIST 



QTR 


11 


12 


13 


14 


15 


16 


17 


18 


19 


20 


S4-1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


84-2 


0 


0 


0 


0 


0 


2 


0 


0 


0 


0 


84-3 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


84-4 


0 


0 


0 


0 


0 


0 


0 


0 


1 


0 


85-1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


85-2 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


85-3 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


85-4 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


86-1 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


86-2 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


86-3 


0 


0 


0 


0 


0 


0 


0 


0 


1 


0 


86-4 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


87-1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


87-2 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


87-3 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


87-4 


0 


0 


0 


1 


0 


0 


0 


0 


0 


1 


88-1 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


88-2 


0 


0 


0 


0 


0 


0 


0 


0 


1 


0 


88-3 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


P(successfu 

repair) 


1.00 


1.00 


1.00 


0.92 


1.00 


0.94 


1.00 


1.00 


0.94 


1.00 
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Table 12. DISPOSAL RATE DISTRIBUTION LIST 



QTR 


21 


A. ^ 


23 


24 


25 


26 


27 


28 


29 


30 


84-1 


0 


0 


0 


0 


1 


0 


0 


0 


0 


0 


84-2 


0 


0 


0 


0 


1 


0 


0 


2 


0 


0 


84-3 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


84-4 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


85-1 


0 


0 


0 


0 


0 


0 


2 


0 


1 


0 


85-2 


0 


0 


0 


3 


0 


0 


0 


1 


0 


0 


85-3 


0 


0 


0 


0 


1 


1 


0 


0 


0 


0 


85-4 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


86-1 


0 


0 


0 


4 


0 


0 


0 


0 


0 


0 


86-2 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


86-3 


0 


0 


0 


0 


0 


0 


1 


0 


1 


0 


86-4 


0 


0 


0 


2 


0 


0 


0 


0 


0 


0 


87-1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


87-2 


2 


0 


0 


0 


1 


0 


0 


1 


0 


0 


87-3 


0 


0 


0 


0 


0 


0 


0 


0 


1 


0 


87-4 


0 


0 


1 


1 


0 


0 


0 


0 


0 


0 


88-1 


0 


0 


0 


1 


0 


1 


1 


1 


0 


0 


88-2 


0 


0 


0 


0 


0 


0 


0 


0 


1 


0 


88-3 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


P(successfii 

repair) 


0.94 


1.00 


0.96 


0.97 


0.97 


0.97 


0.99 


0.99 


0.96 


1.00 
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Table 13. DISPOSAL RATE DISTRIBUTION LIST 



QTR 


31 


32 


33 


34 


35 


36 


37 


38 


39 


40 


84-1 


0 


0 


0 


0 


0 


0 


3 


0 


0 


0 


S4-2 


0 


0 


0 


0 


0 


0 


0 


0 


2 


0 


84-3 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 


84-4 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


85-1 


0 


0 


0 


0 


0 


0 


0 


0 


1 


0 


85-2 


1 


0 


0 


0 


0 


0 


1 


0 


0 


0 


85-3 


0 


0 


1 


0 


0 


0 


0 


0 


0 


1 


85-4 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


86-1 


0 


0 


0 


0 


1 


0 


1 


0 


0 


0 


86-2 


1 


0 


0 


0 


0 


0 


1 


0 


0 


2 


86-3 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


86-4 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


87-1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


87-2 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


87-3 


0 


0 


0 


0 


0 


0 


0 


0 


1 


1 


87-4 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


88-1 


1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


88-2 


0 


0 


0 


0 


0 


0 


0 


0 


0 


- 0 


SS-3 


0 


0 


0 


0 


0 


0 


0 


0 


1 


0 


P(successfu 

repair) 


0.98 


1.00 


0.99 


1.00 


0.99 


1.00 


0.99 


1.00 


0.96 


0.94 
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Table 14. DISPOSAL RATE DISTRIBUTION LIST 



QTR 


41 


42 


43 


44 


45 


46 


47 


48 


49 


50 


84-1 


0 


0 


0 


0 


0 


2 


0 


0 


0 


1 


84-2 


0 


0 


0 


0 


0 


0 


1 


0 


1 


0 


84-3 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


84-4 


0 


0 


0 


1 


0 


0 


0 


1 


0 


0 


85-1 


0 


0 


0 


0 


0 


1 


0 


0 


1 


0 


85-2 


0 


0 


0 


0 


0 


0 


2 


0 


0 


0 


85-3 


0 


1 


0 


0 


0 


0 


0 


1 


0 


0 


85-4 


0 


0 


0 


0 


0 


0 


0 


0 


1 


0 


86-1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


1 


86-2 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


86-3 


0 


0 


0 


2 


0 


0 


0 


0 


0 


0 


86-4 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


87-1 


0 


0 


0 


0 


0 


1 


0 


2 


0 


0 


87-2 


0 


0 


0 


0 


0 


1 


0 


0 


2 


0 


87-3 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


87-4 


0 


0 


0 


0 


0 


0 


1 


0 


0 


1 


88-1 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


88-2 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


88-3 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


P(successfu 

repair) 


1.00 


0.97 


1.00 


0.97 


1.00 


0.99 


0.99 


0.96 


0.92 


0.98 
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B. SAMPLE DATA ANALYSIS 



Table 15. DEMAND DISTRIBUTION ANALYSIS LIST 



No 


Mean 


Mean- 


S.D 


Var 


Q 


1 


0.364 


0.133 


0.674 


0.454 


3.43 


2 


0.545 


0.297 


0.820 


0.672 


2.26 


3 


0.182 


0.033 


0.405 


0.164 


4.97 


4 


1.909 


3.644 


1.700 


2.890 


0.79 


5 


1.909 


3.644 


1.700 


2.890 


0.79 


6 


0.727 


0.529 


1.009 


1.018 


1.92 


7 


0.273 


0.075 


0.647 


0.419 


5.59 


S 


0.455 


0.207 


0.522 


0.275 


1.32 


9 


3.091 


9.554 


2.343 


5.489 


0.57 


10 


1.316 


1.732 


1.293 


2.940 


0.98 


11 


0.895 


0.801 


1.449 


2.099 


2.62 


12 


0.211 


0.044 


0.419 


0.176 


3.99 


13 


0.316 


0.099 


0.946 


0.895 


8.97 


14 


2.632 


6.927 


1.422 


2.022 


0.29 


15 


1.316 


1.732 


1.176 


2.055 


1.19 


16 


0.842 


0.709 


0.834 


0.696 


0.98 


IT 


0.211 


0.045 


0.535 


0.286 


6.36 


18 


1.526 


2.329 


2.118 


4.486 


1.93 


19 


2.579 


6.651 


1.805 


3.258 


0.49 


20 


1.367 


1.869 


1.342 


1.837 


0.98 


21 


1.632 


2.663 


1.862 


3.467 


1.30 


22 


2.421 


5.861 


1.742 


3.035 


0.52 


23 


1.444 


2.085 


1.947 


3.791 


1.82 


24 


2.737 


7.491 


2.642 


6.967 


0.93 


25 


5.362 


31.719 


3.609 


13.025 


0.41 
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Table 16. DEMAND DISTRIBUTION ANALYSIS LIST 



No 


Mean 


Mean^ 


S.D 


Var 


Q 


26 


3.211 


10.311 


2.594 


6.279 


0.65 


27 


45.470 


2067.521 


24.880 


619.014 


0.30 


28 


42.000 


1764.000 


20.700 


428.490 


0.24 


29 


4.684 


21.939 


3.667 


13.447 


0.61 


30 


1.789 


3.201 


1.782 


3.168 


0.99 


31 


8.050 


64.810 


7.920 


62.726 


0.97 


32 


0.421 


0.177 


0.838 


0.702 


5.65 


33 


0.895 


0.801 


0.875 


0.766 


0.96 


34 


5.530 


30.581 


5.230 


27.353 


0.89 


35 


11.160 


124.550 


10.420 


108.576 


0.87 


36 


2.789 


7.779 


2.463 


6.066 


0.78 


37 


1.526 


2.329 


1.467 


2.236 


0.96 


38 


0.947 


0.897 


1.508 


2.274 


2.54 


39 


6.790 


46.104 


5.700 


32.490 


0.70 


40 


3.684 


13.572 


3.284 


10.785 


0.79 


41 


2.053 


4.215 


2.838 


8.054 


1.91 


42 


1.579 


2.493 


2.063 


4.256 


1.71 


43 


0.842 


0.709 


1.167 


1.362 


1.92 


44 


6.053 


36.639 


3.353 


12.496 


0.34 


45 


4.632 


21.455 


3.077 


9.468 


0.44 


46 


5.630 


31.697 


5.380 


30.112 


0.95 


47 


14.890 


221.712 


8.720 


76.038 


0.34 


48 


8.842 


78.181 


3.404 


11.587 


0.15 


49 


3.105 


9.6414 


2.622 


6.8759 


0.71 


50 


9.420 


88.736 


5.370 


28.837 


0.35 
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APPENDIX C. COMPUTER PROGRAMS AND RESULTS 



rrof?ra»a BOKAF ( Inpu t , ou tpu t ) ; 
tyr*3 

tnblO = prrny (1..19.1..0] of roal; 
tab 20 = array (1..01 of roal; 
tobOQ = array [i..0,l..l9,1..0] of real; 
VO r 

clik t boolean; 
l,J,k : Interjcr; 

ln_Lnb : tnbICi; 

DOR^tnli : toblO; 

DHr.Lnb : LnblD; 

nCQ_tob : tab JO; 

OSTO.tnb ; toblO; 

SbQ.tnb : tohlO; 
nO^tnb ; LnblO; 

M.I’_tob : tnbr.O; 

TmiO.tnb : tnbXO; 

FCST^tnb ; tnbVM); 

EXP.tab ; LnblO: 

Din^Lfib : LohlO; 

TEllP_Lob : tobOO; 



PnOCRDUHE tLtlcdi: Integer); 



In 

COSO k of 



1 


11 r 1 to 1 n 


( ' 


.1 t r r 


TUI-lif) AIIAI.YSIS 


); 


2 


u r 1 L 0 1 n 


( * 


1 1 1 ♦ 


nxpoiil-irf I Ab 


SliOUTMIUG 




3 


II r 1 t.n 1 11 


( * 


t t r % 


non 


4 t k k ‘ 


); 


d 


wr 1 Lc 1 n 


( ‘ 


1 M + 


UIIP 


4 ♦ t t t ' 


); 


5 


wr i Lc J n 


( ‘ 


r 1 r t: 


IICQ 


1 4: t > t 


■); 


6 


M r 1 L c 1 II 


( ’ 


» 1 -k 


05TQ 


* t 4- t ' 


■ ); 


7 


w r » 1 f 1 n 


( ' 


•4 t ri; 


SI.Q 


44 ♦ r It 


' ); 


0 


i» r i tc In 


( ' 


t 


no 


t. 4 4- ♦. t. ' 


■ ); 



end (cose U ) 
end ; 



rnOfJKDU UE rcc ( 1 , 3 , k : I n tcgo r ; vo r TEIIP_tnb 
beg ill 

COSO k of 



1 : TUIir.Lnb[k. i . j ] 

2 : Ti:iir.tnli[k. I ,.i ) 

3 : TEIlP_tahtU. i ,J] 

4 : TEMr.tabCk, I .J ] 

5 : TrilP^tobtk, i. j ] 

6 : THllP^tobtU. I , j ] 

7 : TEllP_tnl)(k, 1 . 3 } 

8 : TEllP.tnblk, i . j ] 
end (nose k) 

end : 



;tnb30 ): 



= IH_tnb [i.J]; 

= KXl’^tob [i.j]: 
= DUn.toh [i.j]; 
= UnP^tob [i.j]; 
= nCO_tab [ i » j ) ; 
= OSTO_tnb[ I.j); 
= Sl.n.tab I 1. j J ; 
= UO^tab [I.j); 



rigme 1 1. Conipr.lcr proj^rains for the KAF IModel 
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HAllI PUUGHAll 

begin 

(*>>>> DATA llirt/T IIUUTIIIE <<<<»> 

write In ( ' l)i:ilAlll) DATA !!!*); 
for i := 1 to 19 do 

begin •' 

w r i 1.^ I n ( ’ KOll D ’ » i ) ; 
for J := 1 1.0 0 do 
rend < in_tnb[ i , ,j ) ) ; 
end ; 

wrlteln (' DISPUSAI, RATE !!! '); 

for i ;= I to 19 do 

beg III 

wrlteln ( ’ TOR DIC ,i); 
forj := 1 toOdo 

rend (UlS_tnb [ 
end ; 

(♦>>>> TllEllD AllALYSlS <<<<») 
bog in 

lor i := 1 to 11 do 
beg i n 

TRllU^tnb C I ] : = ( ( ( JH. I.nbf HI . i ) Hll^Lnb( 19 . i J ) ^2)/( lll_tabf IG, i ] > lll_tnb( 17. I Jf 
lH_tnb( 10. i ] 

4 ni..l.nh[ 19.11)): 

if (TRllD..tnh IM >= 0.9) nml (THUD tnl»l i 1 <= 1.1) l.licn 

ALP^tob U j := 0.2 

cine 

ALP.f.nhl ) 1 : = O.'J ; 

end : 

011(1 : 

(♦>>>> EXi-nii i-.iir J Ai. EiHKJTii 1 lu; < < < < n 

for i := 1 toOdo 

KXP.tnb [1.5) := IlL.tnb f l.ij; 
for I := 2 to 19 do 
begin 

for.i ;= 1 toUdo 

EXP_tnb (AI.r..tnb [.}) ^ lH..tab [i.o]) + (( 1 - ALr_Lnb [j)) 

* l•XI’_.tnb ( i-l.j )); 

Olid ; 

for i := 1 to 0 do 

FCST.tob [il := (Al.P.tnb [ij Ill^tob [19. i]) + (<1 - AtP.tnb [i]) 

1 r.Xr.„tnb [I9.i))'; 

(♦>>>> DAILY DF.llAllD RATE <<<c*) 

for I := 1 to 19 do 
bog 111 

forj := 1 toOdo 

DOR^tab t i.j) := lU_tab [ i.j] / 91; 

end ; 



Figure 12. Computer programs for the KAF Model 



86 



(*>>>> 



DEPOT HEPAIR PEnCENT 



<<<<») 



for I := 1 to 19 do 
bc>5 ill 

forj 1 toOdo 

if DlS.tnb [i.j] = U.U then OKP.Lab [i.jj := 1-0 
e l.-e 
b?gin 

if in_tnb Ci,j] = U then DRP_tab [itJ] -= 0.0 
c Is 

UKP^lab fi.j ] := UlS.tab [i.j] / Ul.tab [i.j]; 

**iui : 

if (OlS_lab [i.j] = 1) nncl ( Hl_lab [i.j] = 1) then DHP^tab [i. 
rnd : 



(♦>>>> UCQ K 0:iTq SLO RU <<<<'t) 
fr'r i := 1 to 19 do 
b^r. i II 

for j:= 1 to 0 do 
bnr. i n 

RCQ..tob [i.j] := DOn^tab [i.j] ^ URP^tab [i.j] ^ 120; 

U:'T(}_tal. [i..ij := DOR.tnb [i.j] ♦ ( 1 - DRP.tnb [i.j]) 305; 

Sl.Q^lnb li.jj := r-irt ( 3 ^ (RCO^tab fi.j] + OSTn_tab [i.j])) 
RO^tab [ i . j I := SL<f_tnb[ i . j } i- llCP_tab [ i , j ] + OSTQ_tab [i.jj: 
end : 

eiul ; 



(*>>>> PR HIT lUJUTldK < < < < t ) 

C hk : = fn I !;r» ; 
for k := 1 to n ijo 

b'^f! in 
u r 1 L c 1 n ; 
u r i I e I II ; 

w r i t e 1 II : 

title ( k ) : 

urii.^ Ill ( * cmi 1 2 

H ' ) ; 

wrii.eln ( 

* ) ; 

for i := I to 19 do 
ber; j II 

writc(* 

for j := I ho d tin 

bi'd t '» 

rcc ( i . j . k . TKnp_t nb ) ; 
write (* ’ .TF.HP_tob[ k , i . j ] ) ; 

»'iid ; I 

urihrfn; 
end ; 



Figure 13. Computer programs for the KAF Model 
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j] 0.0: 






ir k 



l Llinn 



bef| in 
w r i tc I n ; 
wr 5 te In ; 

write ( • ALF* ); 

Tor . j ; = 1 to 'I do 

write (* *,ALP_tab [ j ) ) : 
wr i te In ; 
end ; 

i f k - 7. then 

bcf*. in 

w r i le 1 n : 
wr i te In ; 

write ( ' FCT*); 
for j := 1 to A do 

write (' * , FCST_tab[ j ] ) ; 
wr i te ) n ; 
end ; 

w r i Ln 1 II ; 
w r i to I n : 

writoln <* QTM 567 

0 * ) ; 

writcln (’ 

' ) ; 

lor i := 1 to IH do 
bcfl i n 

w r i t. e ( * U * . I ) : 
for .i := 5 to 0 do 
bed i n 

roc ( i . . k . Tl’nr .Lnb ) ; 

wr i l.o( ' ’ ,Tl*.liP,.l.nhl h . i , .] ) ) ; 

end ; 

writcln: 

r.'iid : 

if k = 1 then 

bertiti 

wr i l.e 1 11 ; 
writoln; 
wr itn( * AI.P* ) ; 
r n r .1 : = 5 1. o 0 do 

write<* ’ . AU'^tnb [ ,i J ) ; 
wr j Lein; 
end : 

if k = 7 tlirii 

bed i II 
wr i te In : 
t wr j te J M ; 

writeC FCT'); 
forj :=5to Odo 
write(’ \ FCHT^tnb ( ] ) ; 
writcln: 
emi ; 



end ; 
end . 



Figure 14. Coinpulcr programs for the KAF Model 
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proj[*rom tcrit; 
typo 

nrrny^a = array [)..0J of veal: 
nrrny^K = nrray fl..O) of biLofjer; 
orrny_x = array [1..9J of real; 
var 

nrr_n : orrny^n; 
arr„K : ^nrrny.R; 
nrr.Kc: nrrny.K; 
arr_n : nrray.K; 
nrr.rho : nrray.a; 
nrr.DL: array.ta; 
orr_x : array..x; 
orr.Dn : nrrny.a; 
nrr _Q : orray.n; 
orr^S : array. K; 

5 , totn I . j ,n . i tcm..no, X ,q_tciup : intc^^er; 
fncl . aum.ru . exp,. c:«r;l., pro t_ level ; real; 
do_n;<ajn ; boolean : 

conr;l. 

Biidj'cL 3 lOUdOOOU: *' 

1C = 0.1; 

OA = 0.9: 

(t>>>> STF.P I <<<<:^) 

lunction rtot JU(Kc: integer; rlio: real): real; 

VIII* 

k : iiil.cj^.er: 

\ 

ber.in 

fact ;= 1: 

0: 

for k := 0 to Kc dc) 
bcftiii 

for i:= Kvo doi^iLo Ke-ki 1 do 
fact fact ^ 3/^rho; 

( uriteJiU’k- ‘.I:,* l‘ncL= ’.fact);) 
rnni:r num I (net; 

( wriLcJiK ’rum = *,runi);} 
fan I. : = J ; 
erKi ; 

ncL.l'O := i/auni; 
writcln(‘K) = ’.J/r^um); 
cnrl ; 



Figure 15. 
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procolurc iriit jnl_n( vur nrr_n :nrroy..n) ; 
bcft in 

orr..n[l] ;= 1 JGVf.l. unf 1 . I ; 

Dtr^nf?.] := 3670. nnu. 1 ; 
nrr_j«(3J 5304. UUM. 1; 
nrr^nMl ;= llUOl.Or.tl. l; 
nrr.jil 0) ;= :3;M 13 . (XJ^ 1 . 1; 
nrr_n[Gl := 7307.00^ 1.1; 
nrr_?f7] := llUf;n.72« 1.1; 
nrr_nlOJ := 51043 . 7 I ♦ 1 . I ; 
end ; 



procnlriro clit:piny^rtjnK(Drr_K : arruy.K); 
box ill 



wri i.c )n( ' 


mill 


Kl 


= 


' ,nrr..Kr 1 1 ) 


wr i 1 e ln( * 


PI in 


\:z 


= 


\arr..KLZl) 


wr i In lii( ’ 


mill 


K3 


= 


‘ ,nrr .K13J) 


ur i to ln( ’ 


min 


K4 


r 


*.arr.K(4|) 


wr i tcln( ' 


min 


KO 


- 


■ .nrr^KISJ) 


wr i tc ln( * 


min 


KG 




* »nrr..KlG 1 ) 


wri to Jn( * 


Pi in 


K7 


= 


‘ ,nrr_K( 7 j ) 


wr i te lii( ' 


mill 


KO 


z 


* . nrr_Kl0 j ) 



cnrl ; 



funclioti nl.li power (PL : rnnl;powrr : intr/^er ) : real ; 
bc/t in 

if povmr = 0 tlicii 

iiLli.powcr := 1 
clcn 

nUi .power :r PL y iiLli powcr(0L,powc*r- I ) ; 

rtul : 



runcLi<^ii fnr;(inini : itd.cjlor ) .' iiiLcP.or ; 
bcf? in 

if nnnr - •) Mif»ii 
l :>c := I 
r )rc 

I PC : = mini i I .nr(nnin- J ) ; 

r?iKl; 



rniicLii'ii jict.5(PL: ron I :prol... level : re.nl ) : iiiLcilnr ; 
vnr 

P.r.uni : ten 1 ; 
i : iiiLcp.er: 

: = O.U; 

i':= 0; 

repent < 

r..cnni :r r.niii y MUi_p<iwr»r(UL, i )* (cxp(-I)L))/fac( i ) ; 
i : = i I I ; 

unl.il r.JiUm >= plot.. level; 

:= i*l; 

end ; 



rigure 16 . 
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Figure 17. 



proccOuro (♦cl._ULC vi»r nrr^UL •• orray^nj; 
begin 

orr_l)l.[l) := 0.17/1; 
orr Dl.f71 := 0.032; 
nrr. PI43J := O.KH; 
orr.OM'Jl *•= 0.30*1; 
nrr.DI.l5J := 0.23U; 
orr_DI.[0| := 0.036; 
nrrj)l.l7j :r 0.310; 
orr.ULiOJ := U.GOO; 
end ; 



procrduro dir3ploy..S(ari*_S 
begin 



ur ) tn( * 
wr i te( ‘ 
wri Le( * 
nri Le( ’ 
wr j Lc( ' 
wr i tr( * 
wr i Le( ’ 
wr iteln( ’ 
end ; 



nrrDy_K) ; 



,nrr. 3 l H) 
,nrr_S(2|) 

, arr_ r»( 3] ) 

, orr. 51 'I J ) 
.nrr^S|.3j) 
,nrr_5| 6 1 ) 
.nrr.5|7J); 
\orr_5eni); 



procoJurn gcl_ in J l.( vnr nrr_Kc,.ni i*JJ 
bc(* in 



nrioy.K; vor orr.rlio : arroy.n); 



arr.Kcl 11 


= 


90; 


nrr_Kc[2 1 


r 


1 10; 


or r_Kc( 3) 


r 


103; 


nr r.Kcf '1 1 


r 


103; 


nrr. Kc( 31 


r 


110: 


nrr_Kn( G 1 


r 


100; 


nrr.Kcf?) 


= 


190; 


nrr^ Kr;(.n J 


= 


120; 


nrrJM 1) 


r 


00 ; 


nil* JH 2 1 


r 


no ; 


orr .D( 3 1 


r 


iv:i; 


nrr .1)1 '1 1 


r 


173; 


nrr. I)( 3 j 


- 


00 ; 


nrr_l)| 6 1 


- 


173; 


mrj)f7] 


r 


I/.'); 


nrr.mOJ 


~ 


00 ; 


nrr. rliof 1 1 


= O.UOCGO 


orr..rliol 2] 


1 


= U. 01010 


nrr_r)io[3J 


r O.nuGOl 


orr^rliof'l 1 


= U.OO.T/fJ 


nrr..rl»of 3 


1 


= o.ooorM 


nrr_rlioin j 


= U.0UG02 


orr_rliof7 j 


= n.ni 1329 


nrr_rlio( nj 


= 0.W071 
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nrr_xl 1 J 


- u.uu 


nrr^x[21 


= 0.91 


nrr.x(3] 


= 0.92 


orr_xH ] 


= 0.93 


arr^[5] 


= 0.9'1 


nrr_^[G ] 


= 0.95 


arr^[7] 


= 0.96 


Qrr_x[0 1 


= 0.97 


nrr_^[9 J 


= 0.90 



end; 



bc;|in 

ini tial_a(nrr_m); 

ect^ini t( nrr_Kc , nrr.,I) , oi r_rlio) ; 

repeat 

do.nrtnin := false: 
iLem_no := 0; 
cxp.cost := 0.11; 

repeat 

il.em,.no := item. .no + ); 

ro := get_r0( arr_Kc( ilcmjioJ ,arr_rlio[ item_ncj ) ; 

fnet := 1; 

sum := tl; 

n := 0; 

repeat 

for i:= arr_/%cf i l.eni_iio] clownto nrr_Kc[ i Lem_iio]-n «■ 1 do 
fact: := fact ilem j»o | ; 

fnet := fnet^rO; (nocc-esary to ccmnient) 

write( ' n = * ,n) ; 
sum:= sum + fact; 
writcln( ’ sum = 

fact := 1; 
n := II I- 1; 

uiiU 1(1- sum) <- ( J - (lA) ; 
wr j Leln( *min K = * ,arr. 1)1 i tern, no 1 1 n- I ) ; 
an*. K| i l.eni .nn 1 := nrr.DfilPin nn j + tt -1; 
c.xp.cost := exp.cosL > ar r_n| i tciiijio 1 1 (arr_l)I. itcnijio 1 m-1 ) ; 
until iteni^no = 0; 
if exp_cost <= Ltu'.l/!et. Mien 
d i 5jp 1 ny_ni inK(arr„K) 

• olno 

(lo_nrinin := true; , 

until not do_aj;*aiii; 

(^>>>> SIB’ 11 <<<<:{■) 

wrihelnr prot.levcJ SI S2 S3 SI S3 SC S7 SO'); 

write lii( * ' ) ; 

for X := 1 to U do 
bc(t J n 
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ur i lc( * 


nrr_x(.x J ) ; 


ltcat_no : = 


0; 


gct.DL(arr_ 

repeat 


DU; 


i tcni.no 


: = i Icm^no ^ 1 ; 


t\rr_S[ i tern j»o ] := t»^t_S(nrr J)L[. itcni.no] ,orr_;<(x] ) ; 


until item. 


no - 0; 


d lnploy.S(orr_5) ; 


end; 




(<«>>>> STEl" III <<<<n 


nrr.DoL 1 ] 


= 0.2: 


orr_l)nl.2 j 


= 0.2; 


orr_Un[3] 


= 0.2; 


orr_l)n[^ ] 


= O.'l; 


nrr.D.if 5] 


= 0.2; 


orr. Un[6) 


= 0.2; 


arr_l)a[v j 


= 0.6; 


nr r. I)n| 0 J 
ur iteJn; 


= 0.6; 


wri l.n Jn ; 
ur i telii ( * 


♦ >)>> orriiiAi, n«xiji<iiiii-jrr QUAuriTY <<<<*•); 


for y. := 1 
bcfjin 


to 0 iJo 


nrr.C) lx) := r/irt ((?J (0.1 * (arr..n(.x]/l . 1 ) )t. air.Da [x])/(lC * (arr_atx)/l . 1 ) ) ) : 


write ( ' 


O' .X): 


ur i to ln( 


= ' .nrr_tll y.] ) : 


ri id : 

ur i!.r In ; 
ur i l.e In ; 
ur i I.e Jn ( * 


K 5 « TOTAL): 


ur i l.e In ( * 


•); 


for X 1 


1.0 0 do 


in 




if orr.OIxj <= 2 tlicn Cl temp 
if nrr.CllxJ <= 1 Uien q_temf> : = l: 


tolnJ 


• nir..K(x| ^ arr_..S(xl + q_tcnip; 


ur ile( ’ 


’.X); 


if X = 1 then wrilc( ' ’); 


wril.c( ' 


* .nrr..k| x 1 ) ; 


ur j lc( * 


' ,nrr .SI x | ); 


ur ite( * 


* .q . tnmp) : 


if X = 1 then ur i te ( * ' ) ; 


ur i 1 c( ' 
ur i tcln ; 
t.ikI ; 


■.toi.-il); 


cinJ . 
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Figure 23. Conipulational Results for the KAF Model 
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Figure 24. Computalional Results for llie KAF Model 
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Figure 25. Cominitnlioiial Results for the KAF Model 
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Figure 26 . Compulalional Results for the KAF iModel 
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Figure 27 . Coiiipuintional Results for the KAF Model 
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Figure 29. Comput.ntional Results for the proposed Model (OA = 0.9l)) 
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figure 31. Computational Results for the proposed Model (OA = ().91) 
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Figuio 32. Coinpul.'ilional Results for the proposed Model (OA=0.92) 
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Figure 33. Coiiipu'untiynal Results for the proposed Model (OA= 0.92) 
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